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SECTION  I 

GENERAL  INTRODUCTION 


J 


1.  CCMPOSITE  MATERIALS 


Advanced  composite  materials  have  become  increasingly  Important  in 
recent  years  because  of  their  great  strength  and  relatively  light  weight. 

In  fact,  it  is  estimated  that  the  F-18  aircraft  will  use  advanced  com- 
posites to  the  extent  of  nearly  40%  by  weight  and  80%  outer  surface 
coverage  [1]. 

Advanced  composite  materials  are  laminates,  or  multilayer  "layups", 
of  many  individual  laminae.  A single  lamina  consists  of  a planar  array 

•j* 

of  fibers  (of  boron/boron  tungstate,  graphite  , etc.)  in  an  epoxy  matrix. 
Typical  lamina  thicknesses  are  in  the  range  0.1  - 0.2  mm,  so  that  a 
typical  aircraft  skin  panel  is  made  up  of  at  least  ten  layers.  The 
arrangement  of  the  layers  can  be  varied  to  suJt  the  strength  requirements 
of  the  laminate  in  given  directions.  Typically,  however,  a 0°-90°  or 
0®-45®-90®  layup  will  be  used,  so  that  the  laminate  is  nearly  isotropic 
in  its  mechanical  characteristics  in  directions  parallel  to  its  surfaces. 

A review  of  the  literature  has  shown  that  the  basic  electrical  para- 
meters (e.g.,  permittivity  and  conductivity)  of  advanced  composites  are 
only  now  beginning  to  be  studied  [1,2].  However,  it  seems  generally  to 
be  agreed  that  graphite  composites  behave  as  good  conductors  (in  the  sense 

that  conduction  currents  dominate  displacement  currents  in  the  material) 

4 

with  conductivity  roughly  one-fourth  that  of  pure  graphite  (7.14  x 10 

mho  m ^) , and  that  boron-epoxy  composites  behave  as  mildly  lossy  dielectrics. 


i.e.,  a pyrolyzed  organic  fiber  such  as  polyacrilonitrile. 


In  this  report  we  shall  assume  that  these  materials  are  isotropic. 

It  is  clear  from  consideration  of  the  physical  structure  of  composite 
laminates  that  they  are  probably  anisotropic,  but  for  simplicity  we  shall 
for  the  present  ignore  this  aspect  of  these  materials.  In  a related  study 

[3]  it  has  been  shown  that  the  anisotropic  conductivity  of  graphite  com- 
posites is  not  of  prime  importance  in  characterizing  the  material,  the 
transverse  conductivity  (i.e.,  that  in  the  directions  parallel  to  the 
laminate  surface)  being  the  descriptive  parameter  of  principal  interest. 
Furthermore,  the  anisotropy  of  boron-epoxy  laminates  appears  to  be  weak. 

Since  boron-epoxy  composite  is  such  a poor  conductor,  it  provides 
negligible  shielding  against  penetration  by  electromagnetic  fields.  To 
improve  the  shielding  effectiveness  of  a boron-epoxy  composite  laminate, 
a conducting  screen  may  be  embedded  in  one  of  its  surfaces.  In  this 
study  we  shall  assume  that  the  wires  in  the  screen  are  bonded  at  the 
junctions  and  develop  an  extension  of  the  novj  classic  theory  of  Kontorovich 

[4]  in  order  to  take  into  account  the  effect  of  the  presence  of  a dielectric 
layer  on  the  behavior  of  the  screen. 

The  advanced  composite  materials  we  shall  consider,  therefore,  are 

of  two  types:  graphite  composites,  which  are  modeled  as  homogeneous, 

isotropic,  conducting  materials;  and  screened  boron-epoxy  composites, 

modeled  as  dielectric  layers  with  a bonded  wire-mesh  screen  in  one  surface. 

The  conductivity  of  the  graphite  composites  will  be  taken  to  be  around 
4 -1 

1.5  X 10  mho  m (roughly  one-fourth  the  conductivity  of  pure  graphite) 
and  the  relative  permittivity  of  the  boron-epoxy  composites  will  be 
assumed  to  be  in  the  range  4-5  [3].  Typical  laminate  thicknesses  are 
in  the  range  1-3  mm  and  typical  wire-mesh  screen  parameters  are: 


mesh  size: 


wire  radius: 


20  X 20  to  200  X 200  (per  inch) 
0.05  to  0.15  mm 


wire  conductivity:  1.1  x lo^  mho  m ^ (stainless  steel) 

to  3.7  X 10^  mho  m ^ (aluminum) 


2.  ELECTROMAGNETIC  SHIELDING  CONSIDERATIONS 

We  shall  concentrate  in  this  study  on  the  nuclear  electromagnetic 
pulse  (EMP)  shielding  characteristics  of  graphite  and  screened  boron- 
epoxy  composites.  This  restricts  the  frequency  range  of  interest  to 

g 

f < 10  Hz.  It  should  be  noted  that  characterizations  of  advanced  com- 
posite materials  which  are  valid  in  this  frequency  range  may  not  be 
useful,  say,  for  studying  the  interactions  of  advanced  composites  with 
radar  signals.  The  inhomogeneity  and  anisotropy  of  the  materials,  which 
may  be  neglected  for  EMP  studies,  may  become  critically  important  in 
dealing  with  higher-frequency  interactions. 

The  standard  EMP  waveform  which  we  shall  use  in  this  study  is  [5] 


spectra  of  this  EMP  signal  are  shown  In  Fig.  1.2.  The  function  E(t) 
describes  the  electric  field  time  history  of  an  Incident  electromagnetic 
wave,  which  Is  assumed  to  be  plane.  The  associated  magnetic  field  H(t) 

” E(t)/nQ,  In  which  denotes  the  Intrinsic  Impedance  of  free  space. 

Two  geometrical  forms,  the  Infinite  sheet  and  the  Infinitely  long 
cylindrical  shell,  are  of  principal  Interest  to  us.  We  shall  address 
the  problems  of  electromagnetic  wave  transmission  through  an  Infinite 
sheet  of  advanced  composite  and  of  electromagnetic  wave  penetration 
Into  the  Interior  of  a cylindrical  shell  of  advanced  composite.  Both 
types  of  composites  (graphite  and  screened  boron-epoxy)  will  be  considered 
and  both  frequency  and  time-domain  calculations  of  the  transmitted  or 
penetrated  field  will  be  made,  under  the  assumption  that  the  Incident 
field  Is  as  given  above  In  (1.1)  and  (1.2). 

3.  OVERVIEW  OF  THE  REPORT 

In  Section  II  we  shall  develop  a "boundary  connection  operator" 
which  relates  the  electromagnetic  fields  on  either  side  of  a multilayer 
shield  of  more  or  less  arbitrary  shape.  In  addition  to  Its  utility  In 
non-separable  geometries.  It  simplifies  the  analysis  of  separable 
problems  In  that  the  number  of  regions  to  be  explicitly  considered  Is 
reduced.  This  connection  operator  Is  used  In  the  problem  formulations 
In  Sections  III-VI. 

The  Interaction  of  plane  electromagnetic  waves  with  planar  graphite 
composite  laminates  Is  discussed  In  Section  III  and  the  Interaction  with 
cylindrical  graphite  composite  shells  In  Section  IV.  Both  the  frequency 
and  the  time  domains  are  Investigated  (the  latter  via  Inverse  Laplace 
transformation) . 


Screened  boron-epoxy  laminates  and  their  Interactions  with  plane 
electromagnetic  waves  are  considered  in  sections  V and  VI,  the  planar 
case  In  V and  the  cylindrical  case  In  VI.  Again,  both  frequency  and 
time  domains  are  discussed. 

The  report  Is  summarized  and  the  results  are  discussed  In  Section 
VII.  We  also  suggest  some  extensions  of  the  problems  considered  In 
this  report,  as  well  as  some  logical  next  steps  to  be  taken  In  studies 
dealing  with  electromagnetic  applications  of  advanced  composite  materials. 


SECTION  II 


BOUNDARY  CONNECTION  OPERATORS 

1.  INTRODUCTION 

In  order  to  assess  the  degree  of  electromagnetic  shielding  provided 
by  a layer  of  a given  material  of  specified  shape,  one  must  determine  the 
relationship  between  the  electromagnetic  field  on  one  side  of  the  layer 
and  that  on  the  other.  This  fundamental  mathematical  problem  Is  often 
complicated  by  the  fact  that  the  geometry  of  the  shielding  layer  does 
not  correspond  to  one  In  which  the  vector  wave  equation  Is  separable. 
Therefore,  In  order  to  determine  the  relation  desired,  one  must  either 
develop  exact  solutions  of  such  complexity  and  generality  as  to  be  of 
little  use  for  practical  purposes,  or  take  advantage  of  whatever  special 
circumstances  exist  In  order  to  develop  useful,  albeit  approximate, 
solutions  to  the  problem. 

Fortunately,  such  special  circumstances  occur  In  the  type  of  problem 
In  which  we  are  primarily  Interested.  They  are  that 

(a)  the  thickness  of  the  material  layer  Is  usually  much  smaller 
than  either  of  Its  two  principal  radii  of  curvature, 

(b)  the  wavelength  In  the  material  Is  usually  much  smaller  than 
the  wavelength  outside,  and 

(c)  the  shield  la  usually  made  of  a lossy  material. 

We  conclude  from  the  second  and  third  conditions  above  that  the  electro- 
magnetic behavior  of  a shield  Is  largely  a local  phenomenon.  In  that  the 
fields  at  points  ^ closely 

^^1*  ^2’  ^3  coordinates  of  a system  whose  origin  lies  In  one 

surface  of  the  layer,  the  other  surface  being  at  * d.  and 
are  therefore  coordinates  locally  parallel  to  the  shield  surfaces. 


2 2 1/2 

coupled  when  I + (A?2)  ] » d,  the  thickness  of  the  layer. 

Therefore,  the  characterization  of  the  "transfer  characteristic"  of  the 
shield  layer  requires  consideration  not  of  the  whole  shield,  but  only 
of  local  portions  of  it.  By  virtue  of  the  fl^st  condition  listed  above, 
we  may  consider  these  portions  to  be  planar  and  utilize  a planar  model 
to  determine  the  shield  transfer  characteristic.  Having  developed  this 
characteristic,  we  may  then  apply  the  results  to  geometrical  configurations 
which  are  nonplanar  and  for  which  the  shield  characteristics  may  even 
change  with  position,  provided  that  these  changes  are  sufficiently  gradual. 

The  planar  geometry  which  we  shall  consider  in  this  chapter  is  shown 
in  Fig.  2.1.  The  region  0 £ z £ d^  contains  the  shield  material.  This 
material  may  be  arranged  in  layers  of  different  thicknesses  and  electrical 
properties,  and  the  layers  may  be  separated  by  admittance  sheets.  We 
shall  assume  that  the  individual  layers  and  admittance  sheets  are  homo- 
geneous, linear,  and  isotropic.  In  the  following  sections  of  this  chapter, 
we  shall  develop  an  operator  M (the  boundary  Connection  Operator  or  BCO) 

which  expresses  the  relation  between  the  tangential  components  of  the 

+ 

electromagnetic  fields  on  either  side  of  the  layer  . The  BCO's  for  single 
layers  and  admittance  sheets  are  derived  in  section  B,  and  simplified 
forms  of  these  BCO's  applicable  in  special  circumstances  are  discussed  in 
paragraph  3. 

3.  SINGf.E-LAYER  AND  AMITTANCE  SHEET  OPERATORS 

We  consider  first  the  derivation  of  the  single-layer  operator  M^, 
which  relates  the  tangential  components  of  the  electromagnetic  field  at 
z * + d^  to  those  at  z ■ z^  in  a homogeneous  medium  of  permittivity 

^The  relations  to  be  developed  constitute  an  extension  and  modification 
of  the  ABCD  matrix  approach  to  the  problem  [6]. 

10 


-^Ih 


e and  permeability  The  time  dependence  expCJwt)  is  assumed. 

The  relations  which  we  seek  are  most  conveniently  developed  in  terms 
of  the  Fourier  spectra  of  the  electric  and  magnetic  fields.  Let  us  write 
the  fields  as  follows: 


E(x,y,z) 


HOt.y.zl 


rr 


ECk^.z) 


H(k^,z) 


2- 

e d k. 


(2.1) 


in  which  E and  H denote  the  Fourier  spectra  of  E and  H,  a • k * 0 and 

z t * 

r is  the  position  vector.  E and  H are  expressible  in  terms  of  two  functions 
#(z)  and  y(z)  as 


-jk  xi  ; + fk2  + dlk- 

t z * ue  dz  ^ jme  k ^ dz^k 


(2.2a) 


H = -jk  ^ a ^ + — 4^ — [k^  + — ]ia 

t z dz  jmp^ 


2 2 

^ ® and  y satisfy  the  equation 


(k^  - k )f  - 0 
dz^  ^ ^ 


(2.2b) 


(2.3) 


Solving  eq.  (2.3)  for  $ and  y and  inserting  the  solutions  in  eq.  (2.2) 
yields  expressions  for  the  tangential  field  components  (denoted  by  subscript 


t)  as  functions  of  z'  = z - z^^; 


Et(z') jk^  ^ ®d>  k z*) 


(oT  \ sin  k^z') 


(2.4a) 


If  the  medium  has  finite  conductivity  o,  e should  be  replaced  by  e + ju/a. 
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C2.4b) 


I 

i 


} 


r 

4 

1 


V 

h 

f 

I 

? 


Hj.(z')  = -Jk^  X a^(A^  sin  k^z'  + cos  k^z') 
k 

+ k (A.  cos  k z'  - B.  sin  k z') 
o>p  t ^ z 9 z 


The  constants  A.,  B.,  A,,  and  B,,  which  are  the  arbitrary  constants 
9 9V  V 

associated  with  the  solutions  of  eq.  (2.3),  are  to  be  evaluated  in  terms 
of  the  tangential  field  components  at  z'  * 0,  and 


2 2 - - 
k = k - k -k^ 
z t t 

(2.5) 

We  have,  setting  z'  = 0 in  eq.  (2.4), 

1,(0)  . -Jk,  X ^ k,A,, 

V (2.6a) 

k 

H CO)  = -jk  X i B,  + — k A. 
t t z V (op^  t 9 

(2.6b) 

from  which  we  obtain 

% = ^2  • "t(0> 

k k 
z t 

(2.7a) 

uip 

\ ■ °2  “t  • 

^ k k 
z t 

(2. 7b) 

3k, 

(2.7c) 

3k, 

(2.7d) 

2 - - 

where  k^  B k^*k. . 
t t t 
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FIT 


Now  substitute  eq.  (2.7)  into  eq.  (2.4)  and  set  z'  “ (z  ■ Zj^  + 
We  find,  after  some  simple  manipulations,  the  compact  relation 

fi.cd,) 


1 cos  k d.  U sin  k d. 

Z Z Jw 


-U  sin  k d.  I cos  k d. 

Z 1 Z I 


Ej.(0) 


nii^CO) 


(2.8) 


where  n “ denotes  the  characteristic  impedance  of  the  medium,  I Is 

the  Identity  operator,  and 


U - - a X (k^k^  + k^T) 

jkk  z ^ t t z ' 


(2.9) 


Eq.  (2.8)  is  the  relation  which  lias  been  sought.  We  define  the  single- 
layer connection  operator  as  follows: 


Mj^(Zj^;k;z2) 


so  that  In  general. 


I cos  k^(Zj^  - z^)  U sin  k^(Zj^  - Z2) 


- U sin  k^Czj^  ~ Z2^  ^ ” ^2^ 


(2.10) 


it(^i) 

• • 

nN^(zj^) 

= M^iz^lkiz^)  • 

nSj.(z2) 

(2.11) 


The  notation  chosen  for  the  arguments  of  makes  equations  of  the  form 
of  (2.11)  read  naturally  from  left  to  right.  Note  also  that  M^^(zj^;k;z2)  * 

Mj^(z2;k;zj^). 

We  now  consider  the  connection  between  tangential  field  components 
across  an  Interface.  Such  an  Interface  may  consist  of  a boundary  between 
two  different  media,  or  of  a sheet  Impedance,  or  both.  We  shall  construct 
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'•f*.  9 “k  *-^*^'*'  ^ 

I ■ 1 m * - - ' K 


i 


the  connection  operator  for  the  general  case;  then,  if  the  media  on  either 
side  of  the  boundary  are  identical  or  if  the  admittance  sheet  Is  absent, 
the  connection  operator  will  reduce  to  a simpler  form. 

Let  the  boundary  be  located  at  z = z^.  Then  we  require  that 

lt(2  = z^+)  = E^Cz  * Zg-)  = i^(z  = Zg)  (2.12a) 


lHj.Cz  = Zg  +)  - Hj.(z  = z^-)]  = . i^(z  = z^)  (2.12b) 

is  an  admittance  operator  which  relates  the  surface  current  density 
at  z « Zg  to  the  tangential  electric  field  there: 


•^s  " ^s  • ®t^%>  (2.13) 

It  follows  immediately  from  eq.  (2.12)  that  if  the  characteristic  impedance 

of  the  medium  on  the  "z  +"  side  of  the  sheet  is  q.  and  that  on  the  "z 

» + s 

side  is  n > then 


i,.(z  -) 
t s 

I 

0 

n ^ 

1 

1 

K 

ft 

0> 

1 

a X q Y 
z — s 

— I 

V J 

(2.14) 


Eq.  (2.14)  is  the  relation  which  has  been  sought.  We  define  the  2x2 
inatrlx  in  eq.  (2.14)  to  be  the  admittance-sheet  connection  operator 


M and  write 
s 


a X q Y 

z - s 


— I 


(2.15) 


so  that 


(2.16) 


n.H^CZg-) 


The  notation  chosen  for  the  arguments  of  makes  equations  of  the  form 
of  (2.16)  read  naturally  from  left  to  right.  Note  also  that  M^^(n  “ 

Now  as  an  Illustrative  example,  let  us  construct  the  boundary  connection 
operator  for  the  configuration  shown  In  Fig.  (2.2).  This  Is  a material 
layer  of  parameters  n and  k In  the  region  0 z ;i  d , having  an  admittance 


sheet  on  the  surface  z = 0.  The  medium  outside  the  structure  Is  free 
space,  and  we  wish  to  connect  the  free-space  fields  at  z = 0-  to  the  free- 


space  fields  at  z = d^+.  We  have  the  following  relations: 


E^(O-) 


n^Ht(o-) 


Ej.(0+) 


n H^(0+) 


(2.17a) 


E^(Q+) 


nHj.(Q+) 


= Mj^C0;k;d^)  • 


Et(d,-) 


nHt(d^-) 


(2.17b) 


nHt(d,-) 


= M^(n;0;nQ)  • 


(2.17c) 


from  which  it  is  apparent  that 


in  which 


“t  " ‘ M^CO;k;d^)  • MgCniO;!!^)  (2.19) 

One  will  note  that  the  arguments  in  eq.  (2.19)  read  naturally  from  left 
to  right,  and  correspond  to  the  physical  features  of  the  shield  structure, 
as  is  shown  in  Fig.  2.2.  Furthermore,  note  that 

‘ • MgCn;-Yg;n^)  (2.20) 

Generalizations  to  arbitrarily  complicated  shield  structures  are  obvious. 

In  many  practical  problems,  the  connection  operators  which  have  been 
developed  in  this  section  can  be  simplified  under  certain  conditions.  We 
consider  this  problem  in  the  next  section. 

3.  SPECIAL-CASE  FORMS  OF  THE  CONNECTION  OPERATORS 

In  this  section  we  consider  the  forms  taken  by  the  connection  operators 

when 

(a)  the  electromagnetic  field  is  either  parallel-polarized  or 
perpendicular-polarized 

(b)  the  (effective)  permittivity  of  the  shield  material  is  large 
compared  to  that  of  the  surrounding  medium  (usually  free  space) 

We  also  consider  the  conditions  which  must  hold  for  an  electrically  thin 
layer  to  be  modeled  as  an  equivalent  admittance  sheet. 
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The  tangential  fields  and  H^.  can  be  resolved  into  components 
parallel  and  perpendicular  to  ic^  as  follows:  denoting  either  or 
by  F^,  we  have 

Ft  = (2.21a) 


(2.21b) 

(2.21c) 


in  which  denotes  the  perpendicularly  polarized  part  of  and  denotes 


where 


I cos  k Cz-i  - z-i) 
z 1 2 


X I sin  k^(z^  - Z2) 


M'^(2j^;k;z2)  = 


- ^ X I sin  k^(Zj^  - I cos  “ ^2^ 

z 


(2.23b) 


The  operators  and  are  substantially  simpler  than  and  will  be 
useful  in  dealing  with  shielding  problems  in  which  the  applied  field  can  be 
separated  into  its  parallel-  and  perpendicular-polarized  components. 

We  also  consider  the  operator  in  those  cases  in  which  the  admittance 
operator  diagonalizes  for  perpendicularly-polarized  and  parallel-polarized 
fields: 


— Y’***  I 

s s 


(2.24) 


For  such  cases,  the  connection  operators  M are 


(2.25) 


n Y'  ’"a  XI  — I 
- s z 


When  the  effective  permittivity  of  the  shield  material  is  large 

compared  to  that  of  the  surrounding  medium  (which  is  usually  free  space) , 

the  condition  k^  « |k|  is  valid  for  values  of  k^  corresponding  to  real 

angles  of  incidence.  In  such  cases,  the  approximation  k = k holds  and 

z 

may  be  used  to  simplify  M^,  and  M|^.  We  find  that 


I cos  k(Zj^  - z^) 


M = M'  = MV  = 
I IL  I 


ja  X I sin  k(z^  - z^) 


-ja^  X I sin  k(Zj^  - Z2)  I cos  k(Zj^  - Z2) 


(2.26) 


-T*  • <— ’ 


M 


We  conclude  this  section  by  considering  the  conditions  under  which  an 
electrically  thin  layer  may  be  modeled  as  an  equivalent  sheet  admittance. 
The  connection  operator  for  a layer  In  free  space  Is 


Mf  = M^(n^;0;n)  • Mj^(0;k;d^)  • Mg(n;0;n^) 


(2.27) 


We  define  the  equivalent  sheet  admittance  by  the  relation 


77-1. 


M Cn  ;Y  ;n  ) = M.  (0;k  ;d.)  • 
s o seq  o Jl  ’ o’  I t 


(2.28) 


^1- 


The  factor  C0;k^;d^)  is  Inserted  to  account  for  the  physical  distance 


d^  occupied  by  the  layer.  Expanding  the  product  in  eq.  (2.27)  and  sub- 


stituting in  eq.  (2.28)  yields  the  following  four  equations  expressing  the 
equivalence: 


= = = = 

I = C C.I  + — S S.  U • U„ 
o I n o i o i 


(2.29a) 


0 * c„s  u - — s„c  u„ 

ii  o o n lot. 

O 


(2.29b) 


a X n Y = — c s - s c n 

z o seq  n oil  o £ o 


(2.29c) 


I = C C.l  + s s„u  • u„ 
o I n o £ o £ 
o 


(2.29d) 


in  which 


C = cos  k d„ 
o zo  £ 


C„  = cos  k .d. 
£ zd  £. 


S = sin  k d. 
o zo  £ 


S£  = ‘^zd‘*£ 


(2.30a) 

(2.30b) 

(2.30c) 

(2.30d) 
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(2.30e) 


1 


mmmmm 


U » T, — r — a X (t  it  + k^  I) 
o jk  k z ' t t zo  ' 


o zo 


-1 


Jkk 


zd 


a X (C  k + k^.I) 
z t t zd 


(2.30f) 


and  k^  » k^  - k^,  k^,  = k^  - k^.  Now  If  quantities  of  order  (kd)^  and 
zo  o t zd  t 

higher  are  neglected,  we  find  that  eqs.  (2.29a)  and  (2.29d)  are  satisfied 
as  they  stand;  eq.  (2.29c)  yields  the  well-known  result  [7,8] 


n Y = jk  d,(e  -1)1 
o seq  o i r 

as  we  expect;  but  eq.  (2.29b)  reduces  to 


(2.31) 


° = j Vi 


(2.32) 


k k‘ 

o o 


which  is  not  true  in  general.  It  is  obviously  true  if  k^  = 0,  i.e.,  at 
normal  incidence. 

The  condition  (2.32)  is  also  true  when  the  electromagnetic  field  is 
polarized  perpendicular  to  the  plane  of  incidence.  This  is  so  because 
(cf.  eq.  (2.21c)) 


(iCf^t  ■ ^t^^  * ^t^t  • ° 


(2.33) 


as  is  apparent  by  inspection.  We  now  inquire  under  what  conditions  the 
expression  (2.31)  is  valid  for  parallel-polarized  fields.  We  consider 
the  problem  of  reflection  and  transmission  of  a parallel-polarized  plane 
electromagnetic  wave  by  an  electrically  thin  dielectric  layer  and  by  a 
sheet  admittance  in  free  space.  We  readily  obtain  the  following  results: 

thin  layer:  R"  » [— ^ (e^-l)cos0  - sece(l  - ^)  ] T"  (2.34a) 

o ^r 
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• ^ t Ji- 


-1 


jkd  jk^d 

T"  - [1  + — (e^-l)cos0  + sec0(l  - ^)] 

*^0  r 


(2.34b) 


admittance  sheet:  R"  = — =•  Y"  cos0  T" 
2 s 


(2.34c) 


T"  = (1  + Y"  cos0)"^ 
2 s 


(2.34d) 


2 2 2 

where  k^  = k^sin  0.  It  is  apparent  upon  comparison  of  (2.34a)  and  (2.34b) 
with  (2.34c)  and  (2.34d)  that  If  Y^  Is  given  by  eq.  (2.31),  then  we  must 
require  that  for  the  approximation  of  a thin  layer  by  an  admittance  sheet 
to  be  valid. 


tan  0 <<  |e  I 
' r ' 


(2.35) 


where.  In  general. 


(a  + jtue) 


(2.36) 


We  therefore  conclude  that  the  characterization  of  an  electrically  thin 
layer  as  an  equivalent  sheet  admittance  is  valid  for  perpendicular- 
polarized  fields  generally,  and  for  parallel-polarized  fields  in  which 


k^/k^ 

Is  I 


1 

2 2 

where  k denotes  w u e e ..  being  the  permittivity  of  the  medium 
o o out  out  = r 


(2.37) 


outside  the  sheet  admittance. 
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SECTION  III 

PLANAR  GRAPHITE  C(»1P0SITE  SHIELDS 

1.  INTRODUCTION 

In  this  section  we  shall  consider  the  problem  of  electromagnetic  wave 
transmission  through  a planar  layer  of  graphite  composite  in  the  fre- 
quency and  time  domains.  A brief  description  of  graphite  composite 
materials  has  already  been  given  In  Section  I of  this  report;  we  summarize 
here  certain  facts  pertaining  to  these  materials  which  are  relevant  to 
our  electromagnetic  analysis.  | 

(a)  A layer  of  graphite  composite  Is  a layup  of  anisotropic  laminae  | 

assembled  into  a "cross-ply"  configuration.  Therefore,  the  layer  | 

is  itself  probably  anisotropic  in  behavior.  However,  some  |j 

n 

limited  experimental  data  which  has  been  previously  mentioned  ' j 


i:  ‘ 


i:  i 


h 


; t 


^ I 


[2]  Indicates  that  the  anisotropy  is  weak.  Therefore,  we  shall 
model  the  graphite  composite  as  an  isotropic  material. 

(b)  The  conductivity  of  graphite  composites  appears  to  be  roughly 


one-fourth  that  of  pure  graphite  [1].  This  value  is  7.14  x 10 


-1 


mho  m ; we  shall  use  a value  of  conductivity  equal  to  1.5  x 


4 -1 

10  mho  m for  the  composite  material  in  our  numerical  work. 


(c)  Graphite  composite  is  nonmagnetic  and  we  shall  use  a value  for 
permeability  equal  to  that  of  free  space  in  our  numerical  work. 

In  the  next  section,  we  develop  the  boundary  connection  operator  for 
a planar  graphite  composite  layer  in  free  space.  Then  in  paragraph  3, 
the  problem  of  electromagnetic  wave  reflection  and  transmission  by  such 
a layer  is  solved  in  the  frequency  domain.  The  transmission  of  a transient 
EMP  signal  through  the  layer  is  considered  in  paragraph  4. 


~ \ ' , . - . r-.  -T**  V x-  ^ 

'fi'll- f7rit‘nrf  I f'  ' ..A:.-. -1.  ^ 
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2.  THE  BOUNDARY  CONNECTION  OPERATOR  FOR  A PLANAR  GRAPHITE  COMPOSITE 
LAYER  IN  FREE  SPACE 

The  conductivity  of  a graphite  composite  is  sufficiently  high  that 
the  material  may  be  considered  to  be  a good  conductor  for  frequencies 
throughout  and  beyond  the  EMP  spectrum.  Consequently,  displacement 
currents  may  be  neglected  in  comparison  to  conduction  currents  in  the 
material,  so  that 


2 2 

k = (D  u e - jwu  a -juu  cr 
gg  gg  gg 


(3.1) 


in  which  p , e , and  a denote  respectively  the  permeability,  permittivity, 
g g g 

and  conductivity  of  graphite  composite.  Furthermore,  since  we 


2 2 2 
neglect  in  comparison  to  k in  forming  k^; 

k^d  = kd  = a-j)d/6 


thus 


(3.2) 


where  d is  the  thickness  of  the  layer  and  6 = CZ/uy  o ) 

g g 


1/2 


is  the  skin  depth 

of  the  material.  The  boundary  connection  operator  for  the  graphite  com- 
posite layer  in  free  space  is  thus  given  by  eq.  (2.27),  with  given 
by  eq.  (2.26)  and  k = Cl--j)d/d 


I cosCl-j)d/6 


j X I sin(l-j)d/6 


-j  ^ X I sin(l-j)d/6 


I cosCl-j)d/6 


(3.3) 


in  which  n = j/jmy  /a  is  the  characteristic  impedance  of  the  graphite 
g g 

composite. 

In  the  low-frequency  limit  (d/ 6 « 1) 


-jk  d a XI 
o z 


n.  o d a XI 
o g z 
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(3. A) 


and  when  d/6  >>  1 


f 

i 


Additionally,  we  note  that 


1/2 


(3.5) 

(3.6a 

(3.6b 


3.  PLANE-WAVE  REFLECTION  AND  TRANSMISSION  BY  A PLANAR  GRAPHITE  COMPOSITE 
LAYER 

The  geometry  of  the  problem  to  be  considered  In  this  section  Is  shown 

In  Fig.  (3.1).  A plane  electromagnetic  wave  of  frequency  m Is  Incident 

upon  a planar  graphite  composite  layer  of  thickness  d,  from  the  region 

z < 0;  the  angle  of  Incidence  Is  0.  The  medium  outside  the  layer  Is 

free  space.  The  reflected  and  transmitted  fields  are  to  be  determined 

as  functions  of  the  wave  frequency. 

The  Incident  wave  In  the  region  z < 0 Is  given  by 

-jk  (xslnO  + zcosO) 

^Inc  = ^®0®y  ° (3.7a) 

-jk  (xslnO  + zcos6) 

■^o  ^Inc  “ lEo^a^slne  - i^cos0)  + E|j5y]E  ° (3.7b) 

In  which  Eq  and  Eq  denote  the  amplitudes  of  the  perpendicular-polarized 
and  parallel-polarized  components  of  the  Incident-wave  electric  field. 

The  reflected  wave  In  the  region  z < 0 Is 


Reflected 
Wave  \ 


Transmitted 

'Wave 


Graphite 

Composite 


Incident 

Wave 


Reflection  and  transmission  of  a plane  wave  by  a planar 
graphite  composite  layer:  geometry  of  the  problem 


(3.8a) 


-jk  (xsln6  - zcosO) 

E , - [R'E'a  - R"E"(a  cos0  + a sin0)Je  ° 
ref  Oy  O x z 

-jk  (xsln0  - ZCOS0) 

% “ref  “ lR'EQ(a^cos0  + i^sin0)  + R"EUayjE  ° (3.8b) 

and  the  transmitted  wave  In  the  region  z > d Is 

-jk  Ixsln0  + (z-d)cos0] 

\pa„s  - ^ T’'EJ(VOse  - ;^ai„6)Ja  ° 

(3.9a) 

-jk  [xsln0  + (z-d)cos0] 

n H = [T'E;,(i  sln0  - i cos0)  + T"E"i  ]e  ° 

o trans  O z x Oy 

(3.9b) 

In  which  R',  R",  T',  and  T"  are  the  reflection  and  transmission  coefficients 
for  the  perpendicular-polarized  and  parallel-polarized  parts  of  the  field. 
The  tangential  components  of  the  fields  at  z = 0 and  z = d are 

1^(0)  = (1  + R*)E’  iy 

h!’(0)  = -COS0(1  - r’)e:  i 
O t Q X 

Ej:(0)  = cos0(l  - R")EJJ 
H^(0)  = (1  + R")E|j  iy 

(3.10) 

i^(d)  = T'E^  iy 

Hj;(d)  = -COS0  T'E^  a^ 

E';(d)  = cos0T"E"  I 
t Ox 

"o  ■ ■'"'S  “y 
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Connecting  these  fields  across  the  graphite  composite  layer  using  the 
boundary  connection  operator  given  in  eq.  (3.3)  yields  two  sets  of 


equations  which  may  be  solved  for  the  reflection  and  transmission 


coefficients.  We  have 


-C  - j — COS0  S 1 R'  C - i — cose  S 

^O 

cose  C + j — S COS0  T'  COS0  C - i — S 

n q 


(3.11a) 


'i  1 


cose  c + j — s COS0  R"  cose  c - i — s 

% ^ ^ 

% n 

c + j — COS0  S -1  T"  _c  + j cose  S 


(3.11b) 


in  which  C = cos(l-j)d/6,  S - sinCl-j)d/6.  Solving  eqs.  (3.11),  we  obtain 


expressions  for  R',  T',  R",  and  T"  as  follows: 


R'  js[^cos29 


(3.12a) 


T'  = cose 


(3.12b) 


1 o 

R"  = TJn  jS  — cos  e - — 
D n nj 


T"  = |n-  cose 


(3.12c) 


C3.12d) 


in  which 


fn  9 

D'  = 2cose  C + jS  — cos  0 + — 


(3.13a) 


D"  - 2cos6  C + jS  j — cos  0 + 


Hr. 


2^  n 


(3.13b) 


■ f d-  « * ^ A.  .0  ■ - 1 


Now,  In/n  I = k /n  cr  <<  1,  Neglecting  this  quantity  with  respect 
' o'  o o g 

to  unity,  which  is  consistent  with  neglecting  the  displacement  currents 
with  respect  to  the  conduction  currents  in  the  graphite  composite,  we 
obtain  the  following  expressions  for  the  transmission  coefficients: 


T'  = (C  + 4-  — sec0  S)“^ 

2 n 


(3.14a) 


T"  = (C  + 4-  — cose  S)' 

2 n 


(3.14b) 


Further  approximations  to  the  transmission  coefficients  are  useful.  In 
particular,  it  is  easy  to  show  that  for  the  graphite  composite,  the 
condition  d/ 6 <<  h^c^gd  will  hold  for  frequencies  throughout  and  beyond 
the  EMP  spectrum.  As  a consequence,  the  approximation 


2cose  ,,  d d 


(3.15) 


is  valid  for  all  angles  of  incidence  0,  and  the  approximation 


_ 2sec0  d d 

o g 


(3.16) 


is  valid  for  angles  of  incidence  up  to  a point  only  a few  hundredths  of  a 
degree  from  90°.  Eqs.  (.3.15)  and  (3.16)  show  that  the  effect  on  the 
transmission  coefficients  of  changing  the  wave  frequency  (and  thus  d/6) 
is  almost  completely  separable  from  the  effects  of  changing  ®. 

Therefore,  we  may  express  T'  and  T"  in  terms  of  their  values  at  d/6  = 0 


approximately,  as  follows: 


in  which 


® (1-j)  I csc(l-j)  I 


(3.17) 


« , — .•  V • 


(3.18) 


1 

H 
O- 
1 

2 

nod 

COS0 

= O 
H 

1 

0 g 

sec0 

Numerical  results  pertinent  to  the  analysis  carried  out  in  this 

section  are  presented  in  Figs.  3.2  and  3.3.  In  Fig.  3.2,  and 

Tq  are  plotted  as  functions  of  0 for  various  values  of  the  conductivity- 

4 -1 

thickness  product  a d (note  that  if  a = 1.5  x 10  mho  m , o d = 15d  , 

g g ’ g mm’ 

where  d^  is  the  layer  thickness  in  millimeters) . A plot  of  the  magnitude 

and  phase  of  (l-j)x  csc(l-j)x  as  a function  of  x = d/6  is  given  in  Fig. 

3,3.  Since  the  data  in  Figs.  3.2  and  3.3a  are  presented  in  logarithmic 

form,  the  curves  may  be  used  to  evaluate  T'  and  T"  for  various  combinations 

of  values  of  a d,  0,  and  d/6  simply  by  adding  the  appropriate  quantities, 
g 

It  should  be  noted  that  since  the  magnitudes  of  T'  and  T"  decrease  as  the 
frequency  is  increased,  the  graphite  composite  layer  acts  as  a low-pass 
filter.  This  point  is  important  for  understanding  the  transient  behavior 
of  a transmitted  EMP  signal,  which  is  considered  in  the  next  paragraph. 


4.  EMP  TRANSIENT  PROPAGATION  THROUGH  A PLANAR  GRAPHITE  COMPOSITE  LAYER 
By  using  the  approximate  formulas  for  the  transmission  coefficients 
in  the  frequency  domain  which  were  developed  in  the  previous  section  and 
given  in  eqs.  (3.15)  and  (3.16),  we  may  readily  evaluate  the  temporal 
behavior  of  the  transmitted  wave.  Making  the  formal  substitution  s = jm 
in  T'  and  T",  we  construct  a Laplace  integral  for  the  transmitted  field 
as  follows: 


E,.  (t) 

trans 


■ 2*J  1, 


E„(s) 


nod 
o g 


(3.19) 
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in  which  denotes  the  usual  Bromwich  contour,  x,  is  the  shield  diffusion 
o a 

time,  defined  by 

T , = y a d^  (3.20) 

d g g 

and 

. T = t - ^ [xsin0  + (z-d)cos0]  (3.21) 

c = 1/ the  vacuum  speed  of  light.  In  eq.  (3.19),  "cos0"  is  to  be 
used  when  the  incident  wave  is  polarized  perpendicular  to  the  plane  of 
incidence,  and  "sec0"  is  to  be  used  for  parallel  polarization.  E^(s) 
denotes  the  Laplace  transform  of  the  EMP  signal  discussed  in  Section  i; 


it  is 


in  which 


(3.22) 

(3.23a) 

(3.23b) 

(3.23c)  i 


Furthermore,  it  is  easy  to  show  that  r„  may  be  closed  at  infinity  in  the 

D 

right  half  plane  for  x < 0,  yielding  a null  result  for  E(x),  in  accordance 
with  causality;  and  for  x > 0,  may  be  closed  in  the  left  half  plane  at 
infinity  to  yield  for  E(t): 


ipiinpifP>PV 


I ■ 


^ i 


I 

* I 


Lfe-  . I ... 


E(t)  = 


2A 


n a d 
o g 


fCosO-j 

sec0 


ju(T)  CSC  /ax^  e 


-3t 


- /3t  , CSC  , e 

a d 


» (-1)"  (mr)^(aT  - Bx  ) -n^Tr^x/x  •. 

+ 27  ^ e 

^■2  2 2 2 ® i 

n=l  (n  TT  - aXj)(n  TT^  - 3t^)  ^ 

a a 


(3.25) 


U(x)  is  the  unit  step  function.  The  result  in  eq.  (3.25)  is  based  upon  the 
assumption  that  all  the  poles  of  the  integrand  in  eq.  (3.19)  are  simple, 
i.e.,  that 


(ax^  or  Bx^)  ^ n^Tr^ 


(3.26) 


for  any  n 0.  If  this  condition  is  violated,  the  appropriate  limit  of 

eq.  (3.25)  can  be  taken  to  yield  the  correct  result. 

Plots  of  E(x)  vs.  X,  normalized  by  (2/r)  a d),  are  given  for  the  case 

0 = 0°  in  Fig.  3. A,  for  three  values  of  x^.  The  values  chosen  correspond 

if  is  chosen  equal  to  the  free-space  value  Ait  x 10  ^ henry  m ^ and 
A -1 

o = 1.5  X 10  mho  m , to  thicknesses  d of  1,  2,  and  3 mm.  A non- 

s 

normalized  curve  of  the  incident  EMP  signal  is  included  for  comparison. 
One  will  note  the  increasing  attenuation  of  the  peak  value  of  the  trans- 
mitted signal,  the  increasing  time  delay  before  the  main  buildup  of  the 
transmitted  signal,  and  the  increasing  pulse  width  of  the  transmitted 
signal  as  d is  increased.  All  these  features  are  in  keeping  with  the 
low-pass  filter  characteristic  of  the  graphite  composite  layer  which  has 
been  mentioned  previously. 
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SECTION  IV 


CYLINDRICAL  GRAPHITE  COMPOSITE  SHIELDS 

1.  INTRODUCTION 

In  this  chapter  we  shall  consider  the  problem  of  electromagnetic 
field  penetration  into  the  interior  of  an  infinitely  long  hollow  cylindri- 
cal shell  of  graphite  composite.  A similar  problem  has  been  recently 
addressed  by  Schieber  [9-12],  who  considered  only  cylindrical  shells  which 
were  electrically  thin  (i.e.,  the  thickness  of  the  shell  was  much  less 
than  the  skin  depth  of  the  shield  material) . He  found  that  the  axial 
electric  field  at  the  center  of  the  cylinder  exhibited  strong  peaks  as 
the  frequency  of  the  incident  wave  was  varied,  these  peaks  occurring  at 
the  transverse  resonance  frequencies  of  the  cylinder.  At  these  resonant 
peaks,  the  axial  field  strength  was  equal  to  that  of  the  incident  wave, 
so  that  the  shielding  effect  of  the  cylinder  was  nil. 

The  shielding  behavior  of  a thick  cylindrical  conducting  shell  was 

studied  by  Wu  and  Tsai  [5,13,14].  They  found  that  when  the  cylindrical 

shell  thickness  was  at  least  equal  to  the  skin  depth  of  the  shield 

material  at  the  resonant  frequencies  of  the  cylinder,  the  resonant  peaks 

in  the  axial  field  were  significantly  reduced.  Resonance  effects, 

however,  were  still  evident.  To  assess  the  importance  of  these  resonances, 

Wu  and  Tsai  considered  the  temporal  behavior  of  the  axial  electric  field 

when  the  cylindrical  shell  was  illuminated  by  a normally  incident  TM  - 

z 

polarized  plane  wave  whose  time  dependence  was  that  of  the  standard  EMP 
signal  which  has  already  been  discussed.  They  found  that  the  temporal 
behavior  of  the  axial  field  was  only  slightly  affected  by  the  cylindrical 
resonances  in  a "worst  case"  analysis.  In  their  example,  the  radius  of 
the  cylinder  was  10  m and  the  product  of  the  conductivity  and  shell 


J 


thickness  was  0.1  mho.  The  reason  for  this  behavior  Is,  of  course,  that 
for  cylinders  of  this  size  the  lowest-order  resonances  occur  for  fre- 
quencies far  out  In  the  "tall"  of  the  EMP  spectrum,  and  so  are  only 
weakly  excited.  As  the  radius  of  the  cylinder  Is  decreased,  the  resonant 
frequencies  are  Increased  and  their  excitation  becomes  even  weaker. 

The  cylindrical  shells  to  be  considered  in  this  section  have  radii 
of  2 m or  less  (typical  aircraft-fuselage  sizes)  and  conductivity- 
thickness  products  of  the  order  of  30  mhos.  We  therefore  expect  the 
resonance  effects  on  the  transient  behavior  of  the  internal  fields 
Induced  by  an  EMP  signal  to  be  negligible.  We  introduce  a new  frequency- 
domain  measure  of  shielding  effectiveness  for  hollow  bodies  of  resonant 
size  or  less  which  takes  the  entire  internal  field  into  account.  This 
measure,  based  on  a ratio  of  stored  energies  in  the  interior  of  the 
body,  represents  the  space  averages  of  the  internal  fields. 

In  the  next  paragraph,  the  problem  of  plane-wave  scattering  from  a 
cylindrical  shell  of  graphite  composite  is  formulated  in  the  frequency 
domain,  using  the  BCO  developed  in  the  previous  chapter  for  this  material. 
We  present  results  illustrating  the  energy  shielding  ratio  as  a function 
of  frequency  for  several  representative  cases.  Some  time-domain  con- 
siderations are  discussed  in  paragraph  3 of  this  section. 


2.  FORMULATION  OF  THE  PROBLEM:  FREQUENCY-DOMAIN  RESULTS 


The  geometry  of  the  problem  is  shown  in  Fig.  4.1.  An  infinitely 
long  cylindrical  shell  of  graphite  composite,  whose  inner  and  outer 
radii  are  a and  b,  respectively,  is  illuminated  by  a plane  electromagnetic 
wave.  The  shell  thickness  d = b-a  is  small  in  comparison  to  the  mean 
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Incident  Wave 


Graphite  Composite  Shell 


Plane-wave  penetration  of,  and  scattering  by,  a cylindrical 
shell  of  graphite  composite:  geometry  of  the  problem 


Figure  4.1 


radius  p = The  medium  outside  and  inside  the  shell  is  free  space, 

o 

For  simplicity  we  consider  only  the  case  in  which  the  propagation  vector 
of  the  incident  wave  has  no  axial  component  and  we  shall  Investigate  the 
two  possible  wave  polarizations  separately.  The  total  fields  may  then 
be  obtained  by  superposition. 

Consider  first  the  case  in  which  the  incident  wave  is  polarized  TM 
with  respect  to  z.  The  field  components  present  are  E^,  and 
where 


-ja»p^  y 

(4.1a) 

ay 

(4.1b) 

■ 3P 

1 ^ 

(4.1c) 

P 

i|_  L |ij  ..1  A + A . 0 

P 3P  ( 3pi  .2  - .2  o 


and 

3(|> 

in  the  free-space  regions.  In  the  interior  and  exterior  regions  of  the 
problem,  appropriate  expressions  for  y are 

-E 


(4.2) 


0 < p < a:  y = ^ I A J (k 
— ^ iu)p  ^ n n o 


jup  ' n'"o'‘ 

o n=-" 


(4.3a) 


(4.3b) 


where  A and  B are  to  be  determined,  E is  the  amplitude  of  the  Incident 
n n o 

electric  field,  and 

(b’  = d>  + Tr/2  (4.4) 

o 

is  the  angle  of  Incidence  shown  in  Fig.  4.1. 


1 


In  order  to  evaluate  the  coefficients  and  the  tangential  field 

components  are  connected  across  the  graphite  composite  shell  using  the 

appropriate  connection  operator,  which  is  that  given  in  eq.  (3.3),  modified 

for  the  change  in  the  coordinate  system  (a  is  replaced  by  a ) ; 

z p 


”t  = 


I cos(l-j)d/6 


-j  ^ ip  X I sin(l-j)d/6 
o 


o- 


j -jPp  X I sin(l-j)d/6 


(4.5) 


I cos(l-j)d/6 

The  spectra  of  the  tangential  field  components  at  p=a  and  p=b  are 


E (p=a)  = E A J (k  a)e  a 

t onn  o z 


E^(p.b)  - E^IJ„(k„b)  + '• 


(2)  /I,  “ 


(4.6a) 

(4.6b) 

(4.6c) 


»,(P=b)  - -jEJj;(k^b)  k B„H<^>'(k^b)l.-J"*'  ; 


(4.6d) 

in  which  the  primes  (')  denote  differentiation  with  respect  to  the  argument. 
For  this  problem,  the  spectral  variable  ic^  = (n/p)a^.  Connecting  these 
spectral  components  using  as  given  in  eq.  (4.5)  yields  a pair  of 


equations  for  A and  B : 

n n 


J (k  a)C  + — J' (k  a)S 
no  n n o 

o 


J'(k  a)C  - — J (k  a)S 
n o n n o 


-H^^^(k  b) 
n o 


' (kb) 


k 

J (k  b)' 

n 

no 

B 

J'(k  b) 

n 

no 

(4.7) 


Solving  eq.  (.4.7)  for  A and  B , we  obtain 


1- 


B = ^ (c(J  (k  a)J'(k  b)  - J'(k  a)J^(k  b)] 
n Dl  n o n o no  no 

n '• 

1 

+ S[^  J’(k  a)J'(k  b)  + -^  J (k  a)J  (k  b)] 
n n o no  n no  no  j 


(4.8b) 


where 


D = ClJ’(k  a)Hf^\k  b)  - J (k  a)H^^^’(k  b)] 
n non  o non  o 


- S J ’ (k  ' (k  b)  + ^ (k  a)Hy (k  b)  ] 

nn  on  o nnon  o 
o 


n. 


.(2), 


(4.9) 


C and  S denote  cos(l-j)d/6  and  sin(l-j)d/6  respectively.  and  may 

2 

be  simplified  if  we  assume  that  k^d  « 1 and  |n/n^l  <<  1-  These  condi- 
tions hold  throughout  the  EMP  spectrum  in  the  cases  we  consider;  so  we 
obtain 


nod 


A = [C  + (TTk  p )J  (k  p )H^^^(k  p ) ..  ^ 

n 2 oonoonoo  (l-j)d/o' 


B 


nod  , 

-K  -V-  (Tfk  P„)J:(k  p„) 


(4.10a) 


(4.10b) 


y VES-K  / 7^  • \ , / f. 

o o n'  o o (l“j)d/o 
A useful  frequency-domain  measure  of  the  shielding  effectiveness  of 
the  graphite  shell  is  the  ratio  of  the  time-average  stored  energy  per  unit 
length  inside  the  shell  (0<p<a)  to  the  time-average  stored  energy  per  unit 
length  in  the  same  region  with  the  shell  removed.  Since  this  ratio  is 
obtained  by  integrating  the  stored  energy  density  throughout  the  cylindrical 
volume,  it  represents  an  "averaged"  shielding  effect,  and  so  may  be  more 
useful  than,  say,  the  axial  field  strength,  as  a shielding  effectiveness 
measure  for  a closed  surface.  It  is  easy  to  show  that  the  time-average 
stored  energy  per  unit  length  inside  the  cylinder  is 
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(4.11) 


o o 


In  the  absence  of  the  cylinder,  the  time-average  stored  energy  in  the 


2„2 


same  volume  is  simply  e ira  E /2;  so  performing  the  integrations  in  eq. 


(4.11)  and  normalizing  the  result  yields  the  "energy  shielding  ratio"  r^: 


2 ^ 


= y lA  r [J'(k  a)"^  + J (k  a)"^  1 - 
w ' n'  o ''  no  I ^2^2) 

o 


n=- 


+ 7^-  J (k  a)J'  (k  a)] 
ka  n o n o ^ 
o 


(4.12) 


12  2 

When  ka<<l,  r ==-:r(lA  |^  + I A,  I ),  and  we  obtain  an  approximate  low- 
o w2'o‘  '1‘ 


frequency  expression  for  r^  as  follows: 


w 


4{i 


S . 1-2 

v VI  r\.  .r  , V J / jf  x»n  k p 1 
o g o o (l-3)d/o  o o' 


C - jn  0 d k p 


k a«l 
o 


nod 

+ IC  + j k^p. 


-2, 


2 ■"o^'o  (l-j)d/6'  ) 

Curves  of  r^  as  a function  of  frequency  are  given  in  Figs.  4. 2-4. 4 for 


(4.13) 


various  values  of  p , d,  and  a . The  values  of  p chosen  correspond  to 
o g o 


typical  aircraft;  d = 1,  2,  and  3 mm;  and  a varies  around  the  base  value 

g 


1.5  X 10^  mho  m 


The  analysis  for  the  case  in  which  the  incident  wave  is  polarized  TE 


to  z is  similarly  carried  out.  The  field  components  present  are  H , E , 

z A 


and  E , where 
P 


= jme^  $ 


E = - 
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3p 


(4.14a) 

(4.14b) 


E 

0 o 9(1) 


(4.14c) 


and 


i 

P 3p 


1_  ^ 
p2  3^2 


+ = 0 
o 


(A. 15) 


in  the  free-space  regions.  $ is  given  in  the  interior  and  exterior  regions 
by 


E ® 

0 < p < a;  * = I A'J  (k.  p)e^*'^*"'^' ^ (4.16a) 

— ik^nn  o 

J o n=-“ 

b<p<co;  i = ^ I [J^(k  p)  + p)]eJ*'^*"'*’*^  (4.16b) 

Jo  n=-» 

in  which  and  are  to  be  determined.  Applying  the  connection  operator 
given  in  eq.  (4.5)  to  the  field  representations  obtained  from  eqs. 
(4.14)-(4.16) , we  solve  for  A^  and  B^,  obtaining 


A' 

n 


(4.17a) 


B'  = ^ •(clJ  (k  a)J'(k  b)  - J'(k  a)J  (k  b)] 
n D’  ( n o n o none 

n '• 

+ S[—  J'(k  a)J'(k  b)  + ^ J (k  a)J  (k  b)]!  (4.17b) 

n n o n'o  n n o n'o  j * 

o ' 

where 

D'  = CU'Ck  a)H^^^(k  b)  - J (k  a)H^^^'Ck  b)] 
n n o n o n o n '^o 

- sA  J'(k  a)H^^^'(k  b)  + ^ J (k  a)H^^^(k  b)]  (4.18) 

rinon  o nnon  o 

o 

2 

Approximate  expressions  for  A^  and  B^,  valid  if  k^d  « 1 and  (h/n^l  « 1 
as  we  assume,  are 
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-- 


(4.19a) 
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o o n o o n 


Vo-'  (l-j)d/6 


-1 


('  » -A'  (’fk  D^)J' (k  p^)' 


The  energy  shielding  ratio  for  this  polarization,  r^,  is  given  by 

r’  = ? |a'1^  IJ'(k  a)^  + J (k  a)^(l  - -^1 

w ^ ' n'  ‘ n o no  ( v2  2) 


n=-® 


k a 
o 


(4.19b) 


+ T J (k  a)J  (k  a)  ] 

k a n o n'  o ' 
o 


(4.20) 


1 0 *y 

In  the  case  k^a  « 1,  r^  - ^ ^l^o'  ^ l^ll  obtain  an  approximate 

low-frequency  expression  for  r^  as 

n P d ^ -2 


k a<<l 
o 


■ 2 ^ 2^  ^0^0  (l-j)d/6l 


+ IC  - j 


o g 

2k^p^  (l-j)d/6' 


(4.21) 


Curves  of  r^  as  a function  of  frequency  are  given  in  Figs.  4. 5-4. 7 for 

the  same  values  of  a,  d,  and  a used  in  the  curves  for  r calculated 

g w 

previously. 


It  will  be  noted  that  when  d/6<<l,  r -►laskp  ->-0  but  r' 

w o o w 2 

as  k p -♦■0.  This  is  so  because  the  TE  electric  field  (which  is  normal 
o o 

to  the  cylinder  axis)  terminates  on  surface  charges  on  the  cylinder 

surface.  Thus  r'  contains  only  a contribution  from  the  Internal 
w 

magnetic  field,  while  both  electric  and  magnetic  fields  contribute  to 

r . One  will  also  note  from  eqs.  (4.13)  and  (4.21)  and  Figs.  4. 2-4. 7 
w 

that  the  principal  "break  frequency"  for  both  r^  and  r^  is  approximately 
l/irp  pod. 

0^0  g 
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3.  TIME-DOMAIN  SHIELDING 

In  the  frequency  domain  the  ratio  of  the  Internal  magnetic  field 
to  the  Incident  magnetic  field  is  readily  shown  at  low 


frequencies  to  be 
H 


H 


inc 


nod  -1 

= fC  + j 2 *^0^0  (l-j)d/6^ 


(4.22) 


for  either  polarization.  We  now  construct  a Laplace  integral  for  the 
internal  field  as  follows  (the  incident  waveform  is  the  standard  EMP 
double  exponential) : 

St 


u _ 2dA 

int  27Tj 


^ f (J: L.) 

TTj  J „ s+a  s+3 


ds 
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2d 


'ST,  slnh  /sT,  + — cosh  /st , 
d dp  d 

o 


(4.23) 


in  which  all  quantities  have  been  previously  defined.  Now  2d/p  <<  1 in 

o 

the  cases  which  we  consider,  so  that  the  poles  of  the  integrand,  in 

addition  to  those  at  s = -a  and  s = -3,  are  approximately  located  at 
2 2 

s = -2d/p  T and  -n  tt  /t  (n  > 1).  Thus  we  obtain  for  the  internal 
o a a — 

magnetic  field  the  approximate  expression 
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^o'^o  /3t  , sin  /Bt  , /ax  , sin  /ax  . 

a d d d 

2 2,, 

„ -n  IT  t/x, 

(-1)  e 


+ 2(3-a)x^  I 


1 / 2 2,,-  2 2, 
n=l  (ax,  - n TT  )($t,  - n tt  ) 
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1 

ax  j 


-(2d/p  X )t 
. ° '*  1 U(t) 


(4.24) 


in  which  we  have  made  use  of  the  assumption  ax^  >>  2d/p^.  This  assump- 
tion is  easily  satisfied  in  the  cases  we  consider. 
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. -,2. i. ■ ...  .iA-w 


For  late  times,  the  dominant  contribution  to  (t)  is  the  last 

int 


term  in  (4,24)  (cf.  ref,  [15]) 

H.  . (t) 
int  n p 


o o 


^ L_L-(2d/p^T^)t 


(4.25) 
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4 -1 

mmifo  =1.5x10  mho  m and  u 
8 8 


The  radius  p is  taken  to  be 
o 


A plot  of  ^^^^^^(t)  vs.  t for  three  values  of  is  shown  in  Fig.  4.8 
The  three  values  of  chosen  correspond  to  thicknesses  d of  1,  2,  and  3 

u 

( 

1 m in  each  case.  The  contribution  of  the  expression  in  (4.25)  at  late 
times  is  readily  apparent. 

Comparing  Fig.  4.8  with  Fig.  3.4,  one  will  observe  the  "integrating" 
or  "low-pass"  character  of  the  cylindrical  geometry  Itself.  Thus  an 
incident  signal  is,  in  a sense,  doubly  low-pass  filtered  in  passing  into 
the  interior  of  a cylindrical  region  shielded  by  a conducting  shell. 
Contrasting  behavior  is  evidenced  by  screened  composite  structures,  which 
we  shall  consider  in  the  next  two  sections. 


The  reader  is  reminded  that  in  Cal  Comp  plots  such  as  Fig.  4.8,  the 
axis  labeling  is  such  that  the  exponential  factor  (in  this  case  lO'^) 
is  to  be  associated  with  the  ordinate  values,  not  with  the  axis  label. 
Thus  the  maximum  ordinate  in  Fig.  4.8  is  0.0025. 
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SECTION  V 


PLANAR  WIRE-MESH  COMPOSITE  SHIELDS 


1.  INTRODUCTION 

The  characterization  of  a wire-mesh  composite  layer  which  we  shall 
develop  in  this  chapter  is  that  of  a dielectric  layer  with  an  admittance 
sheet  on  one  of  its  surfaces.  The  presence  of  the  dielectric  layer  in 
the  vicinity  of  the  wire-mesh  screen  is  necessary  for  the  correct  charac- 
terization of  the  screen,  since  its  periodic  structure  causes  the  screen 
to  carry  with  it  a reactive  field  which  "hugs"  the  screen  and  decays 
exponentially  in  directions  normal  to  the  screen.  Thus,  a part  of  this 
reactive  field  penetrates  the  dielectric  layer  and  is  affected  by  its 
presence.  The  equivalent  sheet  admittance  operator  which  describes 
the  screen  is  therefore  dependent  upon  both  the  dielectric  properties 
and  the  thickness  of  the  composite  layer  of  which  it  is  a part. 

The  behavior  of  wire  mesh  screens  in  free  space  has  been  studied  by 
Kontorovlch  [4]  and  others  [16-18]  using  the  so-called  "method  of  averaged 
boundary  conditions."  This  method  is  known  to  produce  good  results  for 
screens  made  of  fine  wire  and  having  spatial  periods  much  smaller  than 
the  wavelength.  More  recently.  Hill  and  Wait  [19]  have  considered  wire 
mesh  screens  using  Fourier  expansions  of  the  wire  currents  and  have  found 
a computational  technique  which  enables  one  to  calculate  the  wire  currents 
relatively  rapidly  when  the  mesh  size  is  small  with  respect  to  the  wavelength. 
We  shall  employ  a method  similar  to,  but  an  extension  of,  that  used  by 
Hill  and  Wait  to  speed  the  convergence  of  the  numerical  solution  for 
the  wire  currents. 
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The  screens  which  have  been  used  to  shield  boron-epoxy  composites 
are  not  generally  fine  enough  to  allow  the  unconcerned  use  of  the  thin- 
wlre  approximation  which  is  basic  to  the  present  theoiry  of  wire-mesh 
screens.  We  shall,  therefore,  compare  the  results  of  our  main  analysis 
(which  is  based  upon  the  thin-wlre  approximation)  with  the  results 
obtained  by  considering  the  limit  in  which  the  screen  is  replaced  by 
a conductor  with  a periodic  distribution  of  small  rectangular  apertures 
in  it.  Using  this  comparison  procedure,  we  can  then  establish  bounds 
on  the  behavior  of  "thick-wire"  mesh  screens. 

In  paragraphs  2 and  3 of  this  section,  we  formulate  the  problem  of 
plane-wave  scattering  by  a fine  wire-mesh  screen  in  the  surface  of  a 
dielectric  layer  of  finite  thickness.  An  equivalent  sheet  impedance  for 
the  screen  is  derived  in  paragraph  4,  and  the  fields  transmitted  through 
the  wire-mesh  composite  layer  are  considered  in  paragraph  5.  The  "per- 
forated screen"  model  is  discussed  in  paragraph  6 and  comparisons  drawn 
between  the  two  models. 

2.  FORMULATION  OF  THE  PROBLEM 

The  geometry  of  the  problem  is  shown  in  Fig.  5.1.  A dielectric 
layer  of  relative  permittivity  (which  may  be  complex)  and  relative 
permeability  1 is  located  in  the  region  0 ^ z ^ d.  A rectangular  grid 
of  conducting  wires  is  embedded  in  the  upper  (z=0)  surface  of  the  layer: 
wires  parallel  to  the  y-axis  are  located  at  x=pa  (p  = 0,+l,...)  and  wires 
parallel  to  the  x-axis  are  located  at  y = qb  (q  = 0,+l,...).  The  wire 
junctions  are  assumed  to  be  bonded.  The  wire  radius  r is  assumed  to  be 
small  in  comparison  to  the  spaclngs  a and  b and  to  the  wavelengths  in 
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Figure  5.1.  Reflection  and  transmission  of  a plane  wave  by  a dielectric 
layer  with  embedded  wire-mesh  screen:  geometry  of  the 
problem 


free  space  and  in  the  layer  material;  thus  the  wire  currents  are  purely 
axial  and  the  thin-wire  approximations  are  valid.  The  wire  impedance 
per  unit  length  is  given  by 


n I (t  r) 
7 = w o w 

w 2Trr  I (t  r) 

L W 


(5.1) 


in  which 


(5.2a) 


T = (jmp  a ) 
w w w 


(5.2b) 


and  0^  denote  respectively  the  wire  permeability  and  conductivity,  and 

I (•)  denotes  the  modified  Bessel  function  of  the  first  kind.  The  time 
n 

dependence  is  assumed  to  be  exp(ja)t). 

The  region  outside  the  layer  is  free  space.  A plane  electromagnetic 


wave  is  incident  upon  the  layer  from  the  region  2<0.  The  incident  electric 


field  is  given  by 


-ik  ‘r 
o 


E.  = E e 
inc  o 


in  which  r is  the  position  vector,  and 


k = k (a  sinO  cos())  + a sin0  sin*  + a cos9) 
oo  x y z 


(5.3) 


(5.4) 


E = E . (a  cos6  cos*  + a cos0  sin*  - a sin0) 

OO0X  y z 

+ E^^(-a^  sin#  + a^  cos#)  (5.5) 

E „ and  E , denote  respectively  the  amplitudes  of  the  0-polarized  and 

Ov  0<p 

#-polarlzed  components  of  the  incident  electric  field,  0 and  # are  the 
incidence  angles  in  spherical  coordinates,  and  k^  = 


^ A m i V I *• 


The  wire  currents  I and  1 are  expressed  in  terms  of  space-harmonic 
X y 


series  as  follows: 


= qb)  = e 


-Ik  qb  “>  > -jk  X 
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(5.6a) 


-jk  pa  “ _ -jk  y 

T / \ V T vn-^ 

I„(x  = pa)  = e I I e 
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(5.6b) 


I and  I denote  the  space-harmonic  components  of  the  currents  I and  I 

xm  yn  x y 


respectively,  and 


k * k sinO  cos*  + 2TTm/a 
xm  o r » 


k = k sin0  sin*  + 2Trn/b 
yn  o r f 


(5.7a) 

(5.7b) 


The  electromagnetic  field  is  expressed  in  terms  of  a vector  function  f 
as  follows: 


H = V X y 


(5.8a) 


E = {VV  + k^I]  • ¥ 
jwe 


(5.8b) 


where  f = IF  a + V a and 
XX  y y 


2-  2- 
V ¥ + k.  f = 0 


(5.9) 


2 2 2 2 
with  k = k e and  e = e e for  0 < z < d;  and  k = k , e = e elsewhere, 
or  or  Q 


o o 

The  periodic  nature  of  the  wire  grid  indicates  that  the  various  field 


quantities  will  be  most  conveniently  expressed  in  terms  of  two-dimensional 
space  harmonic  series  as 


CD  00 


F(x,y,z)  = I I F^^(z)-  ™ - y*" 

m— — n=— 


-jk  X -jk  y 
xm 

e e 


(5.10) 


We  thus  obtain  for  the  tangential  (x  and  y)  space-harmonic  components  of 
E and  H the  following  expressions: 


T » • ^-  4’*  • V X-  . » f 

1..  .r:..:  * V;  . 4^-.- 
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(5.11a) 
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in  which 
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= = 2 2 

In  the  dielectric  layer,  P = P.  , k = k e , and  e = e e ; elsewhere 

mn  dmn  or  or 

= = 2 2 

P = P , k = k , and  e = e . The  space-harmonic  components  of  the 
mn  omn  o o 

surface  current  density  due  to  the  currents  in  the  grid  wires  are 


I I 

J 

sum  b X ay 


(5.13) 


We  now  construct  appropriate  solutions  for  y in  each  of  the  three 

mn 

regions  of  the  problem.  We  have 
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0<z<d;  y =B  cos  k , z + C sin  k , z 
mn  mn  zdmn  mn  zdmn 
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(5.14c) 


2 222  2 2 22-- 
where  k = k - k - k and  k.  =ke  -k^  -k  ;A  -D  are 
zomn  o xm  yn  zdmn  or  xm  yn  mn  mn 

vector  functions  of  m and  n to  be  determined;  and  ir  represents  the 

mn 

incident  wave: 


IT  = j(0£  6 6 P ^ « E 

mn  •'omn  omn  ot 


(5.15) 


6,  denotes  the  Kronecker  delta-function  and  E =(I-aa)*E 
k ot  z z o 

[cf.  eq.  (5.5)]. 
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We  are  now  in  a position  to  construct  the  boundary-condition  equations 

from  which  the  unknowns  A - D may  be  determined  in  terms  of  "n  and 

mn  mn  ' mn 

^smn'  ensuring  continuity  of  tangential  E at  z = 0 and  z = d,  continuity 
of  tangential  H at  z = d,  and  the  appropriate  discontinuity  of  tangential 
H at  z = 0,  we  obtain  the  following  system  of  equations; 
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(5.16b) 


— cos  k ,d  P , • B + — sin  k,dPj*C-P  *0  = 0 
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(5.16c) 


k , sin  k jd  R - k , cos  k .d  C - jk  6=0  (5.16d) 

zd  zd  zd  zd  zo 

We  have  dropped  the  (m,n)  subscripts  in  eq.  (5.16)  and  in  the  following, 

for  simplicity  of  notation.  They  will  be  reinserted  where  necessary  for 

clarity. 

It  is  sufficient  for  our  purposes  to  solve  eq.  (5.16)  for  A and  6. 

We  obtain,  after  some  tedious  but  straightforward  manipulations. 
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in  which 
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(5.17b) 


(5.18a) 


(5.18b) 
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(5.22d) 
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a = 1 + 5 (5.22e) 

Y zo 

b = ^ (e  - 1)  (5.22f) 

Y 

In  order  to  determine  the  still  unknovm  wire  currents  and  thereby 
complete  the  solution  to  the  problem,  we  apply  the  boundary  condition  at 
the  wire  surfaces  that  the  component  of  E parallel  to  a wire  axis  is 
equal  to  the  wire  current  times  its  impedance  per  unit  length.  Since  the 
wires  have  been  assumed  to  be  thin,  it  suffices  to  apply  this  condition 
at  the  tops  of  the  wires,  i.e., 

(q  = 0,+l,...)  (5.23a) 

(p  = 0,+l,...)  (5.23b) 


The  conditions  (5.23)  lead  directly  to  the  following  coupled  set  of 
equations: 


Ey(pa,y,-r)  = Z^Iy(x  = pa) 


ooo 


ooo 


ooo 


Equations  (5.24)  can  now  be  solved  for  the  wire  currents,  and  the  solution 
to  the  problem  is  formally  complete.  We  address  the  problem  of  solving 
eq.  (5.24)  in  the  next  paragraph. 


3.  GRID  CURRENTS 


Equations  (5.24)  may  be  rewritten  as  follows 
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Equations  (5.24)  can  be  solved  numerically  for  I and  I ; however, 

xm  yn 

Hill  and  Walt  [19]  have  shown  that  in  the  case  where  the  dielectric  layer 
is  absent  (l.e.,  the  wire  grid  is  in  free  space),  convergence  of  these 


equations  can  be  extremely  slow.  They  found  that  by  incorporating  a 
Klrchhof f 's-law  junction  condition  into  the  solution,  the  convergence 


■ T *.  m 4 ^ ^ 


was  radically  Improved.  Anticipating  a similar  difficulty  in  the  present 
situation,  we  shall  modify  eq.  (5.24)  by  Incorporating  junction  conditions 
which  should  speed  the  solution. 

Consider  the  functions^ 


f (x)  = ^ 

2iTj  n 


(5.30a) 


r,(x) 


-a  v'  1 _-2jn7rx/a 

0 I 9 a 


(Zir)  n=-“>  n 

which  are  shown  in  Fig.  5.2.  f (x)  has  a unit  jump  discontinuity  at 
x=0,  +a,...;  g (x)  has  a unit  discontinuity  in  slope  at  x=0,  +a,...  . 

di 

US  write  for  I (x,y  = qb)  and  I (x  = pa,y); 

X y 

-jk  qb  -jk  X 

I (x,y  = qb)  = e e [Af  (x)  + s g (x) 

X a X a 


(5.30b) 


Let 


+ I i>  ^-2jm7vx/a 
xm 

m=*-“> 


(5.31a) 


-jk  pa  -jk  y 
X xo*^  ■’  yo-' 

Iy(x  = pa.y)  = e e ■' 


[-Af^(y)  + Syg^^(y) 


+ I r e 

^ vn 


-2iTjny/b 


n=-<» 


yn 


] 


(5.31b) 


The  parameters  A,  s^,  and  s^  are  still  to  be  determined;  it  will  be  noted, 
however,  that  Kirchhoff's  current  law  is  automatically  satisfied  at  the 
wire  junctions.  The  expressions  given  in  eq.  (5.31)  are  equivalent  to 


I * I' 


as 


xm  xm  |_2TTjm  ^2,rm)^J 


IQ 


(5.32a) 


^The  primes  on  the  summation  signs  in  eq.  (5.30)  and  elsewhere  indicate 
that  the  term  n=0  is  deleted. 
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(5.32b) 


A 

I = I'  + «-7--  ^(1-6)  (5.321 

yn  yn  |^2Trjn  (2irn)2j  " 

It  is  now  required  that  the  charge  density  on  each  wire  be  continuous 
at  the  junctions,  i.e., 


31 

31 

X 

X 

3x 

" ax 

x*pa- 

31 

-X 

= 

3y  , 

y=qb- 

3y 

x=pa+ 


(5.33a) 


(5.33b) 


'y=qb-  ^ 'y=ql^ 

for  all  p,q.  The  conditions  of  (5.33)  are  obtained  from  the  continuity 

equation.  Imposing  these  conditions  yields  expressions  for  s and  s in 

X y 


terms  of  A; 


I 


If  these  expressions  are  introduced  into  eq.  (5.32),  we  obtain 


The  final  condition  to  be  Imposed  is  that  the  charge  density  on  the 


’x-wires"  be  equal  to  that  on  the  "y-wlres"  at  the  junctions.  This 


requirement  leads  to  the  condition 


which  yields  the  equation  for  A 


no  potential  difierence  between  the  wires  at  a junction.  Kontorovich  [4] 
used  a more  general  condition  than  (5.36)  to  include  the  effect  of  a 
finite  junction  impedance.  Inclusion  of  this  effect  in  our  analysis  can 
be  shown  to  modify  eq.  (5.37)  by  an  additive  term  t,  in  the  first  square 
brackets;  and 

4 = jwKZ^  (5.38) 

where  is  the  junction  impedance  and  K is  a proportionality  constant 
given  by 

K = [q^(x=pa,y=qb)  - q (x=pa,y=qb)]  (5.39) 

pq  ^ 


is  the  junction  potential  difference  and  q^  and  q^  denote  respectively 
the  linear  charge  densities  on  the  "x-wires"  and  the  "ywires."  It  is 
tempting  to  proceed  with  incorporating  this  junction  condition  in  our 
present  analysis.  However,  it  is  recognized  that  by  Initially  formulating 
our  problem  with  the  x-wires  and  y-wires  coplanar,  the  solution  for  the 
currents  will  converge  to  the  result  for  which  C=0,  no  matter  what  value 
of  ^ is  initially  assigned.  Therefore,  we  shall  consider  only  the  bonded- 
junction  limit  in  the  present  work,  and  reserve  the  case  of  finite 
junction  impedance  for  further  study. 

Incorporating  eqs.  (5.35)  into  eq.  (5.27)  and  simplifying,  we  obtain 
a new  set  of  coupled  equations  for  the  "modified"  space  harmonics  I" 


xm 


and  I'  as  follows:  together  with  eq.  (5.37),  we  have 
yn 
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(2,2)  -jk  r 

r „ , , a ■'  zoon 

n . ,,  -2TTr/a^-l  on  e 

W = Z - -^-1 Jln(l  - e ) T— — — 

n w 2Trjme^  ^J^^zoon 


1 y'  [ (2.2)  e ^-2|m|Trr/a 

■ niLc  2jk^^^a  4lT;a¥ 


(5.42b) 


!h«  = _ ^HL_  tod  - 
^ n 27ri(jje 

n=-oo  o 


- j k r n 

, o\  zomn  s o I I „ /u 

y m ^-2lni.r/b  ^5^^2c) 

nLco  2njk^^^^a  4R^ 


y'  ^EE  = _ Je — tod  - 

“ m 2irja)e 


r " j k r ^ T 

-.  /-o  1 \ zomn  t ol  I / 

I h ■ rr^  ® ^ (5.42d) 

icoe  ^ ran  2mjk  b 4 m it 
-j  o m=-“>  zomn 


in  which 


) 

I 

liin 

2k 

xm 

t 

Inl^ 

mn  0 

1 + e 

r 

f 

Xim 

c 

CM 

1 

r = 
n 

[m  ]-►” 

a — K — 

mn  o 

1 + e 

r 

(5.43a) 


(5.43b) 


. . (1.2)  b xm 

s = lim  a — = —p , 

m II  mn  na  a(e  + 1) 
n -x»  r 


(5.43c) 


t = lim  = (5.43d) 

n I I _ mn  mb  b(e  + 1) 

(m|->«>  r 

The  above  results  hold  when  d > 0.  (If  d = 0,  the  correct  results  are 
obtained  by  setting  = 1). 


•r  ' V ^ « •••-.*. -I 


We  are  now  in  a position  to  calculate  the  wire  currents  in  a 


relatively  efficient  manner,  and  from  the  currents  to  construct  all 
the  field  quantities  of  interest  in  the  problem.  Our  principal  Interest, 


however,  is  in  developing  a relatively  simple  characterization  of  the 
screen  embedded  in  the  surface  of  the  composite,  in  the  form  of  an 
equivalent  sheet  Impedance.  This  problem  is  addressed  in  the  next 
paragraph. 


4.  EQUIVALENT  SHEET  IMPEDANCE  OF  THE  GRID 


In  this  paragraph  we  derive  an  equivalent  sheet  impedance  for  the 


wire  grid  embedded  in  the  surface  of  the  dielectric  layer.  The  relation 


which  we  seek  is 


=Z*J  =Y^*J 
too  s soo  s soo 


(5.44) 


which  relates  the  space-averaged  tangential  electric  field  at  the  grid 
to  the  space-averaged  surface  current  density  [cf.  eq.  (2.16)]. 

From  eqs.  (5.11a),  (5.14a),  and  (5.17a)  we  obtain 


E (z=0)  = -r^ — P • (A  + ir  ) 
too 


— (P  ‘F  ‘A  +P  )*'ir 

JIDE^  OOO  OO  OO  OOO  OO 


-:r  sec0  P • F • r • J 
Z OOO  OO  OO  soo 


(5.45) 


and  using  eqs.  (5.15)  and  (5.25), 


E (z=0)  = b • E :r  sec6  a *5 

too  ot  2 OO  soo 

2k 

o 


(5.46) 
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Now  substituting  eq.  (5.44)  in  (5.46),  we  have 

= = = - '’o  = f 

7 • J = b • E - — K sec0  a • 

S SOO  Ot  -,2  OO  800 


(5.47) 


The  system  of  eqs.  (5.27),  or  (5.37)  and  (5.40),  can,  in  principle  at  least, 
be  reduced  to  a system  of  the  form 


or,  equivalently. 


UI  + VI  = R 
xo  yo  o 


XI  + WI  = S 
xo  yo  o 


bU 

aV  I /b 

xo 

bX 

aW  I /a 

yo 

J U 

(5.48a) 

(5.48b) 


(5.49) 


= b • E 


(5.50) 


and  defining 


t 

-V  • 

bU 

aV 

Z = 

8 

. 1 

bX 

aW^ 

we  obtain 


(5.51) 


Z • J = b • E ^ 
g SOO  ot 

Substituting  eq.  (5.52)  into  eq.  (5.47)  yields 

= = "^o  = 

Z = Z sec0  a 

s g ,.2  OO 


(5.52) 


(5.53) 


• . ■ ' T i»  « I A.  ■ " V • % 

'wM . v#r  * -jr  » ^ i . Jl.  * ^ 


which  is  the  relation  required.  It  is  evident  that  the  problem  of  finding 


the  equivalent  transfer  Impedance  has  been  reduced  to  determining  the 

impedance  Z . This  can  always  be  done  numerically,  of  course:  calculating 
S 

the  currents  I and  I from  (5.37)  and  (5.40)  for  given  E „ and  E ^ 
permits  the  evaluation  of  U,  V,  W,  and  X;  however,  we  recognize  that 
in  the  applications  we  have  Intended  to  consider,  the  mesh  dimensions 
will  always  be  much  smaller  than  the  wavelengths  involved.  As  a 
consequence,  an  approximate  analytic  solution  can  be  obtained. 

The  approximate  solution  is  based  on  the  fact  that  in  the  low- 
frequency  limit,  the  system  of  equations  (5.40)  together  with  the 
auxiliary  condition  (5.37)  has  proven  to  converge  using  only  one  modified 
space-harmonic  in  and  1^  (discussion  of  the  appropriate  numerical 
results  is  deferred  to  paragraph  7 of  this  section).  Therefore,  in 
this  limit,  we  have 


U = U + k f P 
o xo  o 


V = V - k f P 
oo  yo  o 


X = X - k f Q 
oo  xo  o 


W = W + k f Q 
o yo  ^o 


in  which 


f = 


(1  + 


1,2  ,2  .2  2 
kb  , k a , 

_Z£_)  + i (1  + -y  -)i-l 

12  ^ a 12  ^ ^ 


(5.55) 


Z = 
s 


sxx 

sxy 

Z 

syx 

syy 
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(5.56)' 


Thus , if 


then 


, ( 

j 

\ 

« » 


-jk  r 
■'  zoon 


sxx 


co^  - zoon  ( ,, 

° ^ Ijk  {%n’ 

-*  o n=-‘»  ■’  zoon  '■ 


1) 


- k 


xo 


b f rJ^  + 

^ ^iTjn  / 


-fk 


I 


.(1.2)  e 


sxy  jtije  ^ on 
o n=-<» 


-jk  r MU 

zoon  , ik  b 

[J^  + .y°  -] 

2jk  ^2njn  ,2^ 

•’  zoon  (2Trn) 


-f  k <»> 

, xo  r 

Z = r > a 

syx  joie  ' mo 


-jk  r 

/ o 1 \ zomo  .. 
(2,1)  e , 1 


o m=-'> 


2jk 


zomo 


jk  a 
-XO  ■] 

^TTjm  .2^ 


(5.57a) 


(5.57b) 


(5.57c) 


(2TTm)' 


oo 

Z = Z a - — I 

syy  w 3a)e  ^ 


-jk  r 
zomo 


2jk 

o m=-“  ■’  zomo 


|a^2. 

( mo 


2) 


_ ^^XO^  , (2.1)] 

■yo  b ^^27rjm  ^ / 


- k ® + 


(5.57d) 


In  obtaining  the  above  results  we  have  used  the  approximation 

exp(jk^r  COS0)  = 1.  In  order  to  simplify  the  results  of  eq.  (5.56) 

further,  we  take  the  limit  as  k a,  k b,  k.a,  k.b,  k a,  and  k b all 

o o d d xo  yo 

approach  zero,  retaining  only  terms  of  first  order  in  the  final  results. 


(5.58d) 


k s k d k 

^syy  ” V * f'o  ^ ha  ' IT  * i^>ha 

o 


In  which 


. I . ^-.ay/,,-1 

a-l  ” 


(5.59a) 


“>  -2mrr/q  1 + t /e 


L,  = I 

nil 


[ 


nq  r 


1 + i t (e  + — ) 
2 nq  r 


] 


(5.59b) 


..  ..1.  /2mrd. 

t = tanh  ( ) 

nq  ' q 


(5.59c) 


where  q = a or  b. 

Let  us  now  consider  the  eigenvalues  of  the  impedance  matrix  for 
the  case  of  a square  mesh,  a=b.  It  is  easy  to  show  that  the  two  eigenvalues 


Z^^  and  Z^2  ^re  given  by 


k a 


hi  ■ V * 


k a k: 

^s2  “ ^sl  ■ ^2a 

o 


(5.60a) 

(5.60b) 


2 2 2 2 2 

where  k^  = k + k = k sin  6.  Consideration  of  the  eigenvalue  equation 
to  xo  yo  o o 1 


(Z  - Z^.I)  • J = 0 
s si  soo 


(5.61) 


for  i = 1,2  reveals  that 


for  Z = Z , : 
s si 


k^  • J =0 
to  soo 


(5.62a) 


for  Z = Z 

s s2 


k X J * 0 
to  soo 


(5.62b) 


where  k a ^ -j.  ^ _ ^Tie  conditions  of  eq.  (5.62)  are  immediately 

to  xo  X yo  y -i  v ^ 
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f I I . P,  ! 


It  # 


• 1 


r' 


f 

-t 


I 


i 


recognized  as  those  corresponding  to  perpendicular  and  parallel  polarization 
respectively;  therefore 

Z = Z'  = 1/Y'  (5.63) 

sis  s 

is  the  equivalent  sheet  impedance  for  a perpendicular-polarized  field, 
and 

Z _ = Z"  = 1/Y"  (5.64) 

sZ  s s 

is  the  equivalent  sheet  impedance  for  a parallel-polarized  field. 

Several  features  of  these  results  are  worthy  of  some  discussion. 

These  are: 

(a)  Z*  and  Z"  (or  Y'  and  Y")  are  isotropic,  i.e.,  they  do  not  depend 

s s s s 

upon  the  angle  (p.  This  result  is  not  surprising,  since  a bonded 
wire  mesh  in  free  space  is  known  to  behave  isotropically  [4,19] 
at  low  frequencies. 

(b)  Z'  and  Z"  are  not  equal  except  when  = 0;  thus  the  wire  mesh 
is  fundamentally  different  from  a simple  surface  admittance 
layer  in  this  respect. 

(c)  Z^'  is  dependent  upon  so  the  wire  mesh  is  spatially 

dispersive.  This  result  also  agrees  with  the  known  results 
for  the  wire  mesh  in  free  space. 

(d)  zy  depends  upon  and  d.  Physically,  this  occurs  because  for 
parallel  polarization  a normal  component  of  electric  field  is 
present  which  penetrates  the  dielectric  layer.  The  real  impor- 
tance of  this  result  is  that  the  behavior  of  the  grid  depends 
upon  its  surroundings.  As  a consequence,  the  total  behavior 

^ of  the  grid  and  the  dielectric  layer  cannot  be  found  by 
considering  these  components  separately. 
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(e)  The  function  is  bounded  as  follows: 


. -2Trr/a.-l  ^ 2 

)ln(l  - e ) > L„  > 


- 2a  - 1 + e 


ln(l  - 


The  upper  bound  follows  from  setting  = 1 and/or  d = 0 in 
(5.59b);  the  lower  bound  follows  from  letting  d/a  become  large 
in  that  equation. 

(f)  With  regard  to  the  observation  made  in  5 above,  it  must  be 

emphasized  that  je^l  cannot  be  too  large,  since  the  condition 
(k^a/2Tr)  <<  1 must  hold  in  order  for  our  approximate  analysis 
to  be  valid. 

The  thickness  of  the  dielectric  layer  can  be  simply  accounted  for 
by  defining  an  equivalent  relative  permittivity  as  follows: 


eqr 


2 , -2Trr/a.-l  , 

: Ji,n(l  - e ) - 1 

"2a 


(5.65) 


will  depend  upon  d/a,  e^,  and  (weakly)  upon  r/a;  and 
k a k^ 

^ 7 „ j.  <’—2 \n  5® iB_/T  2’n'r/a^— 1 


- V 3"o<2r->ii  - 72 


-]£n(l  - e 


k‘(l  + e ) 
o eqr 


(5.66) 


We  have  calculated  curves  of  e vs.  d/a  for 

eqr 


where  1 < e < e . 

— eqr  — r 

various  values  of  and  r/a;  the  results  are  shown  in  Fig.  5.3.  It 

is  evident  from  these  numerical  results  that  when  d/a  is  greater  than 

unity,  - £^.5  and  since  this  is  the  situation  which  normally  prevails 

in  wire-mesh  composite  applications,  we  shall  set  e = e in  the 

eqr  r 

remainder  of  this  report. 

A generalization  of  this  result  is  obvious:  if  a wire  grid  is 
embedded  between  two  dielectric  layers  of  relative  permittivity  and 
^r2  thicknesses  of  both  of  the  layers  exceeds  a,  then  the  equiva- 

lent surface  transfer  impedances  of  the  grid  will  be  simply 


„ , ^ „ -27rr/a,-l 


jww  a 0/1 

oil  o , o „ -27rr/ax-l  ,- 

Z"  = Z a + — £n(l  - e ) [1  - 


to 


^ ^ 43  1 

s w 2ir 


k (e  , + e 
o'  rl  r2 


-] 


(5.67a) 

(5.67b) 


The  physical  reason  behind  these  results  is  simply  that  the  reactive  field 
of  the  wire  grid  decays  exponentially  in  the  directions  normal  to  the 
grid  as  exp(-2ir  | z | /a) ; and  if  the  layer  boundaries  are  sufficiently  far 
removed  from  the  grid,  these  boundaries  are  not  "felt"  by  the  reactive 
field. 

It  is  convenient  to  express  the  equivalent  sheet  Impedance  for  a 
square  mesh  in  dyadic  form  as 


Z^  = Z + Z - 
s si  s2 


(-jk,o)(-jk,„) 


(5.68) 


or  in  dyadic  operator  form  as 


= - _ V V 

Z = Z T + Z - 
s si  s2  , 2 


(5.67) 


in  which 


k a 

Z = Z a + jn  L, 

si  w o 2ti  la 


k a 

•7  _ • O T 

Z = jn,  L„ 


(5.70a) 


(5.70b) 


"sZ  •'''o  4ir  “2a 

and  Vj.  denotes  the  transverse  or  surface  del  operator.  The  eigenoperators 
for  perpendicular  and  parallel  polarization  are  then  given  by 


Z'  = Z , 
s si 


Z"  = Z , + Z ^V^/k^ 
s si  s2  t o 


(5.71a) 

(5.71b) 
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We  have  used  the  fact  that  is  equivalent  in  the  space  domain  to 
in  the  Fourier  spectral  domain  to  construct  the  result  (5.69)  from 
(5.68)  and  (5.71b)  from  (5.67b).  Curves  of  J^g  shown 

plotted  as  functions  of  frequency  for  typical  screen  parameters  in 
Fig.  5.4. 

In  the  next  paragraph  we  shall  consider  the  problem  of  transmission 
of  a plane  electromagnetic  wave  through  a planar  layer  of  the  wire-mesh 
composite,  using  the  equivalent  sheet  impedance  which  has  been  developed 
in  this  paragraph. 

5.  THE  TRANSMITTED  FIELD 

Having  developed  an  equivalent  sheet  impedance  operator  to  characterize 
the  wire-mesh  screen,  we  may  make  use  of  the  boundary-connection  formalism 
developed  in  Section  II  to  solve  the  problem  of  transmission  of  a plane 
electromagnetic  wave  through  a planar  wire-mesh  composite  shield^.  The 
boundary  connection  relation  for  our  problem  is  (cf.  eq.  (2.19)) 


E,(0-) 


Ej(d+) 


«tn  «t22 


(5.72) 


in  which 


= I cos  k jd 
til  zd 


(5.73a) 


— U sin  k ,d 

tl2  ^ — zd 

r 


(5.73b) 


We  could,  of  course,  calculate  the  transmitted  field  directly  from  the 
analysis  in  paragraph  2 of  this  section.  However,  in  the  frequency  range 
of  interest  for  EMP  studies,  the  approach  based  on  the  equivalent  Z is 
much  simpler.  ® 


(5.73c) 


M „1  = cos  k d a X z + /e  U sin  k ,d 
tzl  o zd  z s r zd 


n. 


M = I cos  k d ^ sin  k .d  (a  x z ^)  • U 

zd  zd  z s 


(5.73d) 


and  U is  given  in  eq.  (2.9)  with  k = k and  k = k 

o r z zd 


For  perpendicular  polarization,  the  elements  M'  are,  when  k ,d  « 1, 

tXj  ZQ 


«;ii  = ^ 


(5.74a) 


^tl2  = ^z  " ^ 


(5.74b) 


1 

";2i  ■ 1 « I 

o 


(5.74c) 


’t\  llT 


m;22  - [1  + JV  n„(zp-Ml 


(5.74d) 


and  for  parallel  polarization,  we  have 


M"  = I 
"til  ^ 


(5.75a) 


”tl2 


- = 

Tr~Z — a X I 
k E z 


(5.75b) 


o r 


”t21  ° I 


or  z 


(5.75c) 


k^,d 


^t22  = + J 


k e %(2s> 


(5.75d) 


o r 


We  have  introduced  the  approximation  k ,d  « 1 because  such  a relation 

Zd 


will  hold  for  a composite  panel  of  this  type  over  the  entire  EMP  frequency 


range:  for  example,  if  e = 4.0  and  d = 2mm,  k ,d  < 0.0084  at  f = 10^  Hz. 

r za  — 
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Now  carrying  out  an  analysis  similar  to  that  done  In  Chapter  111 
we  obtain  for  the  transmission  coefficients  T'  and  T"  the  expressions 


r n ^jkd  •\-l-jkd  cose 

T'  = -jl  + (Zp'^  +-f-  (e^-l)]8ece|  e ° (5.76a) 

r n , jk  d 

T"  = -jl  + [^  (Z;-)'^  + (e^-l)]cos0 


+ sln^0(l  - — )[sec0  + n (Z")~^] 

2 e o s 

r 

to  first  order  In  (k^d) . It  Is  clear  from  eq.  (5.76a)  that  for  perpendi- 
cular polarization,  the  equivalent  sheet  admittance  for  the  wlre-mesh 
composite  layer  Is 


I 

I 


-IK  a COSO 
■’  o 


(5.76b) 


jk  d 

Y'  = (Z')"-^  + — ^ (e  -1)  (5.77) 

s,eq  ' s r ' 

which  Is  In  agreement  with  the  result  given  earlier  In  eq.  (2.31).  For 
parallel  polarization  the  "equivalent"  sheet  admittance  Is 


1 

(Z")'  + -2-  (e  -1) 
s n r 


jk  d 2 1 1 

+ — ^ sin  0 sec0(l )[sec0  + n (Z")  ] 

Ho  o'  s' 


(5.78) 


The  equivalent  sheet  admittances  given  above  can  be  substantially 
simplified  If  we  restrict  our  attention  to  realistic  screen  parameters 
and  the  EMP  frequency  range.  For  example,  a screen  with  a = 0.635  mm 
(40  X 40  meshes  to  the  Inch)  made  of  stainless-steel  wire  (conductivity 
1.1  X 10^  mho  m ^)  of  radius  0.05  mm  yields  a value  of  |ti  (Z')  greater 

O 8 
3 

than  3 X 10  . Thus  we  may  neglect  terms  proportional  to  k^d;  and 


1 


and,  except  for  angles  of  Incidence  6 very  close  to  90°, 


Y"  = (Z")‘^ 
s,eq  s 


(5.80) 


It  is  therefore  apparent  that  the  screened  composite  layer  is  essentially 
equivalent  to  the  screen  alone.  Insofar  as  its  behavior  as  a shield  is 
concerned.  It  should  be  noted,  however,  that  the  effect  of  the  dielectric 
layer  does  Indeed  enter  the  results,  in  that  Z"  depends  upon  the  relative 


permittivity  of  the  layer.  We  give  explicit  expressions  for  ^g’gq 


for  plane  waves  below: 


Y'  = [Z  a + — ~ ind  - e Ipl 

s,eq  w ZTT  * 


(5.81a) 


’'s.e,  “ IV  + ^ ‘“'I  - - rrr))'^ 


(5.81b) 


For  fields  other  than  plane  waves,  sin  0 should  be  replaced  by  k‘'  /k  . 

to  o 


The  transmission  coefficients  T'  and  T"  are  given  approximately  by 


T' 


-jk^d  COS0 


2cos0 


2cos0  + n (Z’)  ^ 
o s 


(5.82a) 


-jk  d COS0 

I’"  ~ ° 


T"  = e 


2sec0 


2sec0  + n (Z")"^ 

O S 


(5.82b) 


These  very  simple  expressions  will  quickly  yield  numerical  values  for  the 
frequency-domain  transmitted  field.  Some  representative  curves  of  the 
magnitudes  of  T'  and  T"  as  functions  of  frequency  for  various  angles  of 
incidence  are  shown  in  Figs.  5. 5-5. 8. 

The  problem  of  transmission  of  a transient  BMP  signal  through  the 
wire-mesh  composite  layer  is  of  interest  and  may  be  easily  attacked. 


2 —1 

For  simplicity,  we  shall  assume  that  Z^  = (irr  a^)  , i.e.,  that  the  screen 


wires  have  radii  small  in  comparison  to  the  skin  depth  over  the  frequency 
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in  which  E^(t)  is  given  in  eq.  (1.1).  The  factors  y' /&'  and  y"/S"  are 


given  by 


^ 2^ 


(5.86a) 


n n r a 
o w 


y"  _ 2a 
■pT  = 2 — 


(5.86b) 


IT  n r a 
o w 


Note  that  the  quantity  Ttr  o^/a  plays  the  same  role  in  these  expressions 

as  does  the  quantity  a d in  related  formulas  for  the  graphite  composite 

8 

case.  Curves  of  Ej.(t)  vs.  t at  normal  incidence  are  shown  for  two 
representative  cases  in  Figs.  5.9  and  5.10. 

The  shapes  of  the  E^(t)  curves  of  Figs.  5.9  and  5.10  can  be 
understood  readily  in  terms  of  the  "differentiating"  or  high-pass  charac- 
ter of  the  bonded  wire-mesh  screen.  The  sheet  Impedance  of  the  screen 
comprises  both  resistive  and  inductively  reactive  components,  and  thus 
Increases  in  magnitude  with  increasing  frequency.  The  effect  of  this 
behavior  on  the  transmission  coefficient  in  the  frequency  domain  is 
evident  in  Figs.  5.5  through  5.8:  configurations  with  = 1.1  x 
10^  mho  m ^ yield  values  of  |t' | and  |t"|  which  are  nearly  constant  with 
respect  to  frequency  over  the  range  of  Interest,  and  the  sheet  impedance 
is  primarily  resistive.  Increasing  the  wire  conductivity  to  3.72  x 10^ 
mho  m ^ decreases  the  resistive  component  of  the  sheet  impedance,  making 
the  reactive  component  relatively  more  Important.  For  the  two  mesh 
slze/wlre  radius  combinations  considered,  the  transmission  coefficient 
at  low  frequencies  is  higher  for  the  case  a = 0.212  mm,  r = 0.051  mm  than 
for  the  case  a = 0.635  mm,  r » 0.127  mm,  as  a consequence  of  the  fact 
that  the  second  configuration  contains  more  metal  per  unit  area,  even 


though  the  mesh  size  is  larger.  On  the  other  hand,  the  inductive  reactance 


of  the  0.635/0.127  configuration  is  larger  than  that  of  the  0.212/0.051 
mesh,  for  equal  wire  conductivity;  the  magnitude  of  the  transmission 
coefficient  for  this  case  is  correspondingly  greater  at  higher  frequencies. 

The  effect  of  this  behavior  on  the  time-domain  transmitted  electric 
field  is  evident  in  Figs.  5.9  and  5.10.  The  0.212/0.051  mesh  yields 

-12 

an  E^(t)  response  which  is  nearly  Independent  of  up  to  x * 5 x 10 
sec;  this  early-time  behavior  is  principally  a result  of  the  inductive 
reactance  of  the  sheet  impedance.  At  later  times  the  resistive  component 
of  the  sheet  Impedance  exerts  the  dominating  Influence  on  the  transmitted 
field. 

The  0.635/0.127  mesh,  on  the  other  hand,  yields  a much  higher 

i 

amplitude  for  the  transmitted  field  at  early  times  (again,  up  to  x = 5 x 
“12 

10  sec,  the  response  is  nearly  Independent  of  a ) as  a consequence  of  I 

w i 

its  larger  inductive  reactance.  The  late-time  amplitude  is  again  dominated 

by  the  resistance  of  the  screen,  although  some  undershoot  is  apparent 
-9 

for  X > 2 X 10  sec  for  the  case  of  the  larger  wire  conductivity. 

The  most  Interesting  comparison  to  be  made  among  the  four  curves 

presented  is  that  between  the  two  higher-conductivity  cases.  The  finer 

2 

mesh,  even  though  its  equivalent  conductivity- thickness  product  irr  ff^/a 
is  smaller  by  a factor  of  2 than  that  of  the  larger  mesh,  yields  a 
smaller  transmitted  field.  This  observation  leads  one  to  conclude  that 
the  conductivity-thickness  product  alone  cannot  be  used  as  a measure  of 
the  effectiveness  of  mesh  shields.  The  time-domain  behavior  of  the 
Incident  fields  must  also  be  taken  into  account;  in  fact,  an  interesting 
problem  for  further  study  would  be  to  design  a screen  to  yield  minimum 
transmitted  field  under  a given  set  of  constraints. 
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One  should  also  compare  the  curves  In  Figs.  5.9  and  5.10  with 
those  in  Fig.  5.4  for  the  planar  graphite  composite  layer.  The  pulse 
shapes  have  already  been  discussed.  Of  interest,  however,  is  comparing 
the  peak  transmitted  electric  field  amplitudes  for  the  two  types  of  com- 
posite. For  example,  a 2 mm  thick  graphite  panel  of  conductivity  a = 1.5  x 

4 -1  -5 

10  mho  m yields  a peak  transmitted  electric  field  amplitude  of  9.7  x 10  ; 

the  0.635/0.127  mesh  yields  a peak  value  of  45  x 10~^;  and  the  0.212/0.051  mesh, 

15  X 10  The  screened  panel  would  thus  seem  generally  to  be  a poorer 

shield  than  the  graphite  panel,  at  least  at  early  times.  We  will  return  to 

this  type  of  comparison  in  the  next  section. 


6.  THE  PERFORATED  SCREEN  MODEL  FOR  A WIRE  MESH 

In  this  paragraph  we  briefly  consider  the  problem  of  transmission  of 
a plane  electromagnetic  wave  through  a periodic  square  array  of  square 
apertures  in  a perfectly  conducting  plane,  in  order  to  compare  the  results 
to  those  previously  obtained  for  the  wire  mesh.  One  purpose  is  to  gain 
some  insight  into  the  behavior  of  the  mesh  screen  when  the  assumption 
r/a  <<  1 does  not  apply  and  the  thin-wire  approximation  breaks  down. 

Let  the  screen  be  in  the  surface  z=0  and  let  the  apertures  be  square 
of  side  length  h and  centered  at  (x  = (m  + l/2)a,  y = (n  + l/2)a)  for 
< m,n  < <».  The  geometry  is  shown  in  Fig.  5.11.  It  is  well  known 
(cf.,  e.g.,  [20])  that  an  electrically  small  aperture  in  a conducting 
screen  can  be  modeled  in  terms  of  equivalent  electric  and  magnetic 
dipoles.  We  therefore  consider  the  electromagnetic  field  radiated  by 
a planar  array  of  electric  and  magnetic  dipoles  in  free  space  located 
at  (x  = (nrfl/2)a,  y = (n+l/2)a,  z=0).  The  electric  and  magnetic 
dipole  moments  are  assumed  to  be  given  by 


T * * • 


p (x  = ma  + a/2,  y = na  + a/2)  * p a e 

eoz  z 


-jk  a(iD+l/2) 


-jkyoa(n+l/2) 


(5.87a) 


p (x  = ma  + a/2,  y = na  + a/2)  = p e 

rao 


- j k^^a  (nri-l/  2 ) - j k a (n+1/  2 ) 


(5.87b) 

in  which  p^  denotes  the  electric  dipole  moment  and  p^  denotes  the  magnetic 

dipole  moment,  p and  p are  constants:  furthermore  p = p a + 
eoz  mo  mo  mox  x 


p a . 
moy  y 


When  the  period  of  the  array  a is  very  small  in  comparison  to  the 


free-space  wavelength,  the  radiated  electric  field  is  easily  shown  to  be 


, kp  -ikx-jky-jkz 

„ 1 . xo'^eoz  . xo  yo-^  ■'  zo 

E = — ^ ( P )e 

X „ 2 me  moy 

2a  o 


(5.88a) 


, kp  -ikx-iky-ikz 

E . 1 + p “ yo'  “ (5.88b) 

y „ 2 me  mox 

2a  o 

1 1 o o “jk  x-jk  y-jk  z 

1 r,  - 1 „ 2.,2^,-’xo-^  yo-'  zo 

E = — x [k  p - k p - -= (k  + k )]e 

z 2a^k  xo  moy  yo  mox  ^me^  xo  yo 

zo 


(5.88c) 


2 2 2 2 

in  the  region  z > 0:  also,  k = k - k - k . 

zo  o xo  yo 


Now  the  dipole  moments  p and  p are  given  in  terms  of  the  fields 

•^eoz  '^mo  ® 

E^^  and  which  would  exist  in  the  plane  z-0  if  the  apertures  were 
closed  by  the  relations 


p = jme  a E 
eoz  o e zo 


(5.89a) 


p = -imp  a H 
mo  o m to 


(5.89b) 


in  which  and  denote  respectively  the  electric  and  magnetic  polariz- 
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abilities  of  the  apertures.  Latham  [20]  has  calculated  these  polarlz 


1 + 0.104(h/a) 


1 - 0.185(h/a) 


Let  us  now  assume  that  a perpendicularly  polarized  plane  wave  is 


incident  upon  the  perforated  screen  from  the  region  z<0.  Then  E and 


in  which  E 


denotes  the  electric  field  amplitude  in  the  Incident  wave 


k sln6  sin4i,  and  k 


which  8 and  (J)  are  the  incidence  angles.  The  transmitted  electric  field 


amplitude  is  now  easily  shown  to  be 


When  a parallel-polarized  plane  wave  is  incident  upon  the  screen 


sin9 


denotes  the  electric  field  amplitude  in  the  Incident  wave.  The 


transmitted  electric  field  amplitude  is 


i 

i 

! 

t 

I 

r- 

»• 


t 

••  t 


Examining  eqs.  (5.92)  and  (5.95a)  for  we  see  that  the  expressions 

— Jln(l-e  ^ and  a /a^  (cf.  eq.  (5.90b))  are  to  be  compared.  We 

TV  m 

shall  compare  them  on  the  basis  of  equal  optical  coverage  by  letting 

h = a-2r.  The  two  exp/.essions  are  shown  as  functions  of  r/a  in  Fig. 

5.12.  The  wire-mesh  model  factor  is  always  larger  than  the  perforated- 

screen  factor;  the  former  is  accurate  in  the  limit  r/a  ■*  0,  and  the 

latter  is  accurate  in  the  limit  r/a  ->•  0.5.  It  is  clear  that  the  wire 

mesh  model  overestimates  the  transmitted  field  amplitude,  and  therefore 

underestimates  the  induced  currents  on  the  screen,  for  values  of  r/a 

which  are  not  very  small  in  comparison  to  unity. 

Examining  eqs.  (5.94)  and  (5.95b)  for  we  see  that  the  relevant 

comparison  is  between  the  constant  factor  0.5  for  the  wire-mesh  model  and 

the  ratio  a /a  for  the  perforated-screen  model.  This  ratio  is  shown  as 

e m 

a function  of  r/a  in  Fig.  5.13.  It  is  apparent  that  the  mesh  model 

2 

overestimates  the  importance  of  the  "sin  9"  term  in  the  transmission 
coefficient  for  parallel-polarized  fields.  The  reason  for  this  is  probably 
the  use  of  the  thin-wire  approximation  in  the  analysis  of  the  mesh  model. 

In  effect,  the  equivalent  permittivity  of  the  mesh  for  electric  fields 

normal  to  the  plane  of  the  mesh,  caused  by  Induced  line  dipole  moments 
on  the  wires,  has  been  ignored.  It  is  clear  from  our  earlier  analysis 
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of  the  effect  of  a dielectric  environment  on  the  wire  mesh  that  including 

2 

these  line  dipole  moments  would  decrease  the  coefficient  of  the  "sin  0" 


term  and  thereby  improve  the  accuracy  of  our  approximate  solution  for 


the  sheet  impedance  operator  Z . This  aspect  of  the  problem  would  seem 


worthy  of  further  study 


7. 


NUMERICAL  COMPUTATION  OF  GRID  CURRENTS 


The  equations  (5.37)  and  (5.40)  can  be  solved  for  the  modified 

current  harmonics  I'  and  I'  by  means  of  the  procedure  outlined  in  this 
xm  yn 

paragraph.  We  begin  by  substituting  eq.  (5.37)  into  eq.  (5.40),  thereby 
eliminating  the  variable  A from  further  consideration  and  yielding 


y (U  6 + P k )I'  + y (V  - P k )!'  = R 6 

^ ' m mn  2 m xn  xn  ^ mn  2 m yn  yn  o 


(5.95a) 


y (X  - ^ Q k )I'  + ^ (W  6 + Q k )I'  = S6 

^ nm  2^xnxn  ^ mmn  2^mynyn  oi 


(5.95b) 


as  the  coupled  sets  of  equations  for  I'  and  I'  . Also, 

xn  yn 

1,2  k2  J-  2 

.kb  , k a . 

/I  . yo  X . /I  . xo  ^ T -1 

^ ^b  ^2  24  ^ ^2  24 


(5.96) 


Now  define  matrices  M , M , M , and  M <=uch  that 

XX  xy  yx  yy 

(M  ) = U 6 + ^ P k 

XX  mn  m mn  2 m xn 


(5.97a) 


(M  ) = V - ^ P k 

xy  mn  mn  2 m yn 


(5.97b) 


(M  ) = X - ^ Q k 

yx  mn  nm  2 ^ > 


(5.97c) 


(M  ) = W 6 + ^ Q k 

yy  mn  m mn  2 m yn 


(5.97d) 


and  let  I'  and  I*  denote  the  modified  space  harmonic  current  vectors. 

X y 

Finally,  let  R^  and  denote  the  "driving"  vectors;  these  each  possess  only 
a single  nonzero  element.  Then  the  coupled  equations  (5.95)  may  be  written 


i'  + M • i'  = R 
X xy  y o 


(5.98a) 


I'  + M 

X 


(5.98b) 


' f 


This  system  is  readily  solved  for  and  We  obtain 


where 


!'  * 

N 

• R - 

N • S 

(5.99a) 
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and  the  solution  is  formally  complete.  A can  now  be  calculated  using 
eq.  (5.37). 

The  calculation  of  the  grid  Impedance  Z defined  in  eq.  (5.52)  is 

S 

straightforward.  We  readily  obtain 


(5.101) 

Now  using  eq.  (5.53),  we  may  calculate  the  equivalent  sheet  impedance  Z^. 

The  efficient  numerical  evaluation  of  the  matrix  elements  defined  in 

eq.  (5.97)  has  already  been  discussed  (cf.  eqs.  (5.42)  and  (5.43)). 

We  have  calculated  the  eigenvalues  of  the  equivalent  sheet  impedance 

operator  Z as  a function  of  N,  the  number  of  modified  space  harmonics 
8 

A*  M 

included  in  I*  and  I'  (note  that  N =•  1,3,5,...),  for  the  case  d/a  “ 1, 

X y 


f (N  ) 

b XT  oo 


- (N  ) 
a yr  oo 


h (N  ) 

b xs  oo 


- - (N  ) 
a ys  oo 


-1 
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b/a  1,  r/a  = 0.01,  e = 4.5,  a = «,  0 » 60®,  where  ^ = 0®  and  45® 

r w 

and  k^a  = 1.0  and  0.1.  The  results  of  these  computations  are  shown  In 

Figs.  5.14-5.17,  in  which  values  of  X'  and  X"  (Z  ’ = iX  ’ ) are 

s s s s 

plotted  vs.  N.  The  rapid  convergence  of  X'  and  X"  with  increasing  N 

s s 

is  evident  (note,  however,  that  in  the  (|>  = 45®  cases  this  is  not  at 
all  surprising,  since  there  is  no  discontinuity  in  the  junction  currents 
when  ij)  = 45®,  135®,  225®,  or  315®,  by  symmetry).  The  use  of  the  N=1 
case  to  derive  an  analytical  approximation  for  Z^  when  k^a  « 1 is 
justified  by  the  rapidity  of  convergence  apparent  in  Figs.  5.16  and 
5.17  for  k a = O.l'*'. 


Note  that  for  typical  mesh  dimensions  the  maximum  value  of  k a occurring 
in  the  EMP  frequency  range  is  of  order  10“^. 
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SECTION  VI 


CYLINDRICAL  WIRE-MESH  COMPOSITE  SHIELDS 

1.  INTRODUCTION 

In  this  section  we  shall  consider  the  interaction  of  a plane  elec- 
tromagnetic wave  with  an  infinitely  long  cylindrical  shell  of  wire-mesh 
composite.  As  in  Section  IV,  we  shall  consider  only  the  cases  in  which 
the  propagation  vector  of  the  incident  wave  has  no  component  along  the 
cylinder  axis. 

In  the  next  paragraph,  the  problem  is  formulated  in  the  frequency 
domain  using  the  equivalent  sheet  impedance  developed  in  the  previous 
section  to  describe  the  screen  in  the  surface  of  the  composite  shell. 

Some  representative  frequency-domain  results  are  presented.  In  paragraph 
3,  we  briefly  consider  the  transient  behavior  of  the  internal  magnetic 
field  when  the  circumference  of  the  cylinder  is  small  in  comparison  to 
the  wavelength  over  the  frequency  range  of  interest  and  the  temporal 
behavior  of  the  incident  wave  is  that  of  the  standard  EMP  signal. 

2.  FORMULATION  OF  THE  PROBLEM 

The  geometry  of  the  problem  is  shown  in  Fig.  6.1.  An  infinitely 
long  cylindrical  shell  of  wire-mesh  composite  of  inner  and  outer  radii 
a and  b respectively  is  illuminated  by  a plane  electromagnetic  wave. 

The  medium  inside  and  outside  the  shell  is  free  space  and  the  shell 

thickness  d ■ b - a is  small  in  comparison  to  the  mean  radius 

The  side  of  the  wire-mesh  composite  in  which  the  grid  is  embedded  is 

the  outer  surface  of  the  cylinder.  For  simplicity  we  shall  again  consider 

only  the  case  in  which  the  incident  wave's  propagation  vector  has  no 

axial  component,  and  we  shall  investigate  the  two  possible  wave  polarizations 

separately.  The  total  fields  may  then  be  obtained  by  superposition. 
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Incident  Wave 


Composite  Shell 


Plane  wave  penetration  of,  and  scattering  by,  a cylindrical 
shell  of  boron-epoxy  composite  with  a bonded  wire-mesh  screen 
in  the  outer  surface:  geometry  of  Che  problem 


Figure  6.1 


The  shell  Itself  Is  modeled  as  an  equivalent  sheet  Impedance,  based 


on  the  results  obtained  for  the  planar  case  in  the  previous  section. 


We  shall  denote  the  mean  shell  radius  by  p . Z",  which,  as  will  be 

o s 

recalled,  depends  upon  the  transverse  propagation  constant,  will  be 


denoted  Z"  for  each  angular  eigenfunction  in  the  field  expansion,  and 
s,n 


= il- 

“^to  “-ton  2 

Po 


(6.1) 


Consider  first  the  case  in  which  the  incident  electromagnetic  wave 
is  polarized  TM  with  respect  to  z.  In  this  case,  the  field  components  are 


E = -jmp  y 
z ■'  o 


H = ^ 

■ 3p 


H = i ^ 
P P 3i(> 


(6.2a) 

(6.2b) 

(6.2c) 


where 


-E 


0 < p < p : 
— — o 


P IPq: 


[ A J (k  p)ej^^^-'^’^ 

n n r\  ' 


iwy  ^ n n ' o 


(6.3a) 


-E 


I [J  (k  o)  + B Hf^^(k  p)]e^"^''’”'^'^  (6.3b) 


imp  ^ ■ n'  o 

o n=-® 


n n 


A and  B are  to  be  determined,  E is  the  electric  field  amplitude  of 
n n o 

the  incident  wave,  and  <p'  = (}>  + ti/2,  where  (p  is  the  angle  of  incidence 

o o 

as  shown  in  Fig.  6.1.  The  appropriate  Z for  this  polarization  is  Z'. 

s s 

Now  the  boundary  conditions  at  P'P^^  are  applied.  We  have 


E_(P  * P„-)  = E (p  = p +)  = E (p  = p ) 
z o z o z o 


«A<P  = P«+>  ■ “ (2')“^  E (p  = p ) 

<b  O p O 3 z o 


(6.4a) 

(6.4b) 
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from  which  expressions  for  and  are  readily  obtained.  We  find 

(2)  1-1 

A = 1 + irk  p I J (k  p )H^  ^(k  p ) ^ 

n o o 2 n o o n o o' 


7 

os  , 


J = -irk  p — % J (k  p )A 

n o o 2 n o o n 


(6.5a) 


(6.5b) 


When  the  incident  wave  is  polarized  TE  with  respect  to  z,  the  field 


components  are 


H = jwe  $ 
z o 


3<I> 

E = - — 
<P  3p 


E =ili 

P P 3<t> 


(6.6a) 


(6.6b) 


(6.6c) 


where 


0 < p < p : f = -r^  J A’J  (k  p)e^^^‘*’~'*’'^ 

- - ° iK  n=-co  " " ° 


(6.7a) 


P ^ Po'  ^ " W ^ B^H^^^k^p)]e^"^‘'’"‘^'^  (6.7b) 

o n=-" 

A^  and  B^  are  to  be  determined.  At  p = p^,  the  boundary  conditions  to 
be  applied  are 


= ^*(P  = pQ^ 


(6.8a) 


"e«-  ■ »o*)  - «.<»  ■ »„-)  - ■ Po> 


(6.8b) 


in  which 


%<»„)  ^ I E (»„) 

n=-<»  ’ ^ 


(6.9) 


We  find  for  A'  and  B' 
n n 


" -1 

r n (Z  ) ^ 

At  A , o s,n 

A'  = 1 + Ttk  p ^ J 

n o o 2 


’(k  p )H^^^'(k  p ) 
noon  o o 


(6.10a) 


(6.10b) 
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The  energy  shielding  ratios  r^  and  r^  are  given  in  terms  of  and 
in  eqs.  (4.12)  and  (4.20).  To  obtain  expressions  for  r^  and  r^ 
appropriate  for  the  screened  composite  cylindrical  shell,  A^  and  A^ 
given  In  eqs.  (6.5a)  and  (6.10a)  are  substituted  Into  eqs.  (4.12)  and 
(4.20)  respectively.  When  k a « 1,  approximate  expressions  for  r and 


r are 
w 


w 


k a<<l 
o 


l-jkp  n(Z')^J!,nkp 
o o o s o o 


-2 


n„(Z')‘^,-2^ 


1 + jk  p 
•'  o o 


(6.11a) 
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(6.11b) 


Curves  of  r^^  and  r^  as  functions  of  frequency  are  g’ven  for  various  values 

of  screen  parameters  and  a cylinder  radius  p^  =•  1 m in  Figs.  6.2  and  6.3. 

(a  , rather  than  a,  denotes  the  mesh  size  in  this  section).  As  in  the 
m 

related  curves  for  the  graphite  shell,  we  note  that  unless  a * r -►1 

w w 

and  r^  -»■  0.5  as  f ^ 0.  The  reason  for  this  behavior,  as  before,  is  that 
r^  contains  only  a contribution  from  the  internal  magnetic  field,  while 
both  the  internal  electric  and  magnetic  fields  contribute  to  r^.  Also, 
we  have  not  included  curves  for  r^  in  the  case  = 1;  over  the  frequency 
range  shown,  they  are  indistinguishable  from  the  curves  for  the  case 


4.5. 
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Comparison  of  Figs.  6.2  and  6.3  with  Figs.  4. 2-4. 7 reveals 

that  for  cylindrical  geometries  the  screened  boron-epoxy  composite  shells 

with  = 1.1  X 10^  mho  m ^ are  comparable  to  graphite  composite  shells 

in  shielding  effectiveness,  while  those  with  - 3.72  x 10^  mho  m ^ are 

decidedly  superior.  The  reason  for  this  is  that  the  principal  break 

frequency,  below  which  the  shielding  is  relatively  poor,  is  approximately 
2 2 

l/(iTu  p (nr  o /a)):  and  nr  o /a,  which  is  the  equivalent  conductivity- 

O O W V 

thickness  product  for  the  mesh  is  comparable  to  a d for  graphite  shells 

8 

when  = l-l  * 10^  niho  m but  is  significantly  larger  than  o^d  when 

a = 3.72  X 10^  mho  m 
w 

Not  shown  in  the  r^  and  r^  curves  is  the  higher-frequency  effect  of 
the  screen  inductance.  This  is  apparent  in  the  time-domain  internal  field 
respont>e,  which  is  discussed  in  the  next  paragraph. 


3.  TIME-DOMAIN  SHIELDING 

The  ratio  of  the  internal  magnetic  field  to  the  incident  magnetic 
field  in  the  frequency  domain  is  readily  shown  to  be  given  when  <<  1 


it"  = j Vo  r (2;)  ] 


(6.12 


for  either  polarization.  This  result  does  not  involve  the  relative 
permittivity  of  the  composite  shell  at  all,  since  the  case  in  which 
the  magnetic  field  is  perpendicular  to  the  cylinder  axis  corresponds 
to  perpendicular  polarization  (and  Z'  does  not  involve  e ),  and  the 

S IT 

case  in  which  the  magnetic  field  is  parallel  to  the  cylinder  cixis 
utilizes  only  the  n=0  angular  eigenfunction  in  the  low-frequency  limit 
(and  Z"  _ does  not  involve  e ) . 


When  the  screen  wires  are  perfectly  conducting,  the  result  in 
eq.  (6.12)  reduces  to 

-2Trr/a_ 


“iM  /,  ^ ">o  , -1„ 
inc  m 
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Tip 
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Jln(l  - e 


I . ^ 

nij-l 


(6.13) 
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which  is  a constant,  independent  of  frequency.  Thus  the  time  dependence 
of  the  internal  magnetic  field  is  identical  to  that  of  the  incident 
field,  except  for  a constant  attenuation.  This  behavior  is  attributable 
to  the  combined  "differentiating"  effect  of  the  perfect  screen  and  the 
"integrating"  effect  of  the  cylindrical  geometry. 

Another  simple  result  of  some  Interest  is  obtained  if  we  assume 


that  over  the  frequency  range  of  Interest,  where 


(xr^o  )~^ 
w 


is  the  resistance  per  unit  length  of  the  screen  wires.  In  this  case,  we 


-2Trr/a 


find  that  when  (a  /jrp  )Jln(l-e 
m o 


m.-l  , 

) « 1, 


H.  a L jio  + 2'irR  /p  L 
int  m w o 


Inc 


irp  ju)  + 2R  a /y  p 
o w m o o 
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in  which  L denotes  £n(l-e  ) The  impulse  response  of  the  Internal 

magnetic  field  is  therefore  approximately  given  by 


a L 2R  a •-(2R  a /y  p )t 

H.  ,.(t)  = — 6(t)  + — e w ">  o o 
int'  ’^Po  PoPo 


(6.15) 


so  that  the  Internal  field  response  to  a standard  EMF  waveform  is 
dominated  at  late  times  by  the  slowly  decaying  exponential  function  in 
the  second  term  in  (6.15).  More  specifically,  under  standard  EMP 
excitation,  the  Internal  magnetic  field  is  given  by 
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in  which  we  have  made  use  of  the  fact  that  (B  and  a)  » 2R  a /p  p in 

w m o o 

typical  cases.  Representative  curves  of  ^ shown  in 

Figs.  6.4  and  6.5  for  various  values  of  screen  parameters  and  a 


cylinder  radius,  = 1 m. 

The  curves  in  Figs.  6.4  and  6.5  show  very  clearly  the  enhancement 

in  shielding  effectiveness  with  increasing  wire  conductivity.  In  the 

curves  of  Fig.  6.4,  for  example,  the  peak  value  of  the  internal  field 

drops  from  about  0.09  to  about  0.005  when  a is  increased  from  1.1  10^ 

w 

to  3.72  X 10^  mho  m One  will  also  note  the  long  "tails"  of  the  internal 

field  waveforms  due  to  the  term  in  the  expression  in  (6.16)  depending 

on  t as  exp(-2R  a t/p  p ). 

w m o o 

Comparing  the  lower-conductivity  curves  in  the  two  figures,  one  notes 
immediately  that  the  0.635/0.127  mesh  is  a better  shield  than  the  0.212/0.051 
mesh;  this  is  a consequence  of  the  larger  equivalent  conductivity-thickness 
product  of  the  former  configuration.  It  is  also  apparent  that  the 
"integrating"  effect  of  the  cylindrical  geometry  has  delayed  the  onset  of 
the  peak  signal,  in  comparison  with  the  planar  case  (cf.  Figs.  5.9  and 
5.10).  Thus  the  wire  conductivity  (or,  more  precisely,  the  equivalent 
conductivity-thickness  product)  controls  the  behavior  of  the  Internal 
waveform.  In  the  higher-conductivity  cases  the  inductive  reactance  of 
the  screen  becomes  more  important:  one  will  notice  the  peak  in  the 


132 


.0 

response  for  the  0.635/0.127  mesh  near  t * 2 10  sec.  This  peak  is 

absent  in  the  response  of  the  0.212/0.051  mesh  because  of  Its  lower 
inductive  reactance. 

Finally,  if  one  compares  the  responses  in  Figs.  6.4  and  6.5  with 
those  in  Fig.  4.8  for  the  graphite  composite  cylindrical  shells'  internal 
field,  one  notes  that  the  peak  internal  field  for  the  thickest  graphite 
composite  shell  is  about  the  same  as  that  for  the  "worst-case"  screened 
boron-epoxy  composite  shell;  furthermore,  the  "best-case"  screened  boron- 
epoxy  shell  response  is  lower  than  that  of  the  thickest  graphite  shell  by 
a factor  of  roughly  20.  One  could  conclude  from  these  results  that  a fine 
mesh  screen  made  of  highly  conducting  (A1  or  Cu)  wire  would  be  a much  better 
shield  for  cylindrical  geometries  and  EMP  waveforms  than  even  a thick 


graphite  laminate. 


SECTION  VII 


CONCLUDING  REMARKS 


I 


1.  SUMMARY  AND  DISCUSSION 

In  this  report  we  have  considered  both  the  frequency-domain  and  time- 
domain  shielding  provided  by  graphite-epoxy  and  screened  boron-epoxy  com- 
posite laminates  in  planar  or  cylindrical  geometries,  with  respect  to  EMP 
signals  and  spectra.  Graphite-epoxy  composites  have  been  modeled  as 
isotropic,  homogeneous  conducting  materials  and  boron-epoxy  composites 
as  isotropic,  homogeneous  dielectrics.  The  screen  used  in  screened 
boron-epoxy  panels  has  been  modeled  under  the  assumptions  that  the  wire 
junctions  are  bonded  and  that  the  thln-wlre  approximation  applies.  These 
assumptions  may,  in  fact,  not  be  well-satisfied  in  practice  (the  wire 
junctions  may  be  imperfect  due  to  oxidation  at  the  junction  points  of 
overlapping  wires  and  the  condition  that  the  wire  radius  is  small  in 
comparison  to  the  mesh  period  is  not  generally  well  met);  however,  they 
provide  a reasonable  starting  point  for  analysis  and  yield  conservative 
estimates  of  the  mesh  shielding  effectiveness,  since  the  thln-wlre  approxi- 
mation tends  to  underestimate  the  mesh  current. 

We  have  developed  a boundary  connection  operator  formalism  which  is 
a generalization  of  the  transmission  matrix  of  network  theory.  Such  a 
formalism  will  be  useful  in  dealing  with  shielding  problems  Involving 
bodies  of  non-separable  shape  as  long  as  the  local  radii  of  curvature  are 
small  in  comparison  to  the  shield  thickness.  The  boundary  connection 
matrix  (or  "super-matrix")  possesses  many  interesting  properties  which 
have  not  been  touched  on  here  but  which  will  be  the  subject  of  a forth- 
coming Interaction  Note.  The  principal  utilities  of  the  boundary  connection 
operator  are  its  generality  (in  that  general  multilayer  shields  with  sheet 
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Inmlttances  between  layers  are  easily  described)  and  the  fact  that  Its 
use  effectively  reduces  the  number  of  regions  which  must  be  considered 
In  a shielding  calculation  by  at  least  one. 

The  transmission  of  electromagnetic  plane  waves  through  planar 
graphite  composite  laminates  Is  straightforward  to  analyze,  and  some  useful 
approximate  formulas  for  the  transmission  coefficients  of  such  a struc- 
ture have  been  given.  The  transmitted  time-domain  field  behaves  In  a 
manner  easily  understood  on  the  basis  of  the  low-pass  filter  characteristic 
of  conducting  materials.  We  have  calculated  some  time-domain  curves  for 

i 

t3rplcal  cases  and  have  found,  for  example,  that  the  peak  value  of  the  j 

i 

electric  field  transmitted  through  a 2 mm  thick  graphite  composite  panel 
-4 

is  approximately  10  of  the  peak  value  of  the  Incident  field. 

Penetration  of  electromagnetic  fields  into  the  interior  of  cylindrical 
shells  of  graphite  composite  has  been  considered.  Over  the  frequency  range 
which  has  been  of  principal  Interest  to  us,  and  for  cylinder  radii  on  the 
order  of  2 m or  less,  resonance  effects  are  of  little  importance.  Con- 
sequently, low-frequency  approximations  are  useful.  We  have  suggested  a 
new  frequency-domain  measure  of  the  shielding  effectiveness  of  closed 
surfaces  In  terms  of  the  time-average  electromagnetic  energy  stored  In  the 
Interior  of  the  region,  and  have  given  simple  formulas  for  this  quantity 
for  cylindrical  bodies.  We  have  also  presented  representative  curves  to 
Illustrate  the  dependence  of  this  shielding  parameter  on  frequency  for 
several  typical  graphite-composite  cylinders,  and  have  pointed  out  that 
the  principal  "break  frequency"  below  which  the  shielding  is  poor  is  given 
approximately  by  l/ny^o^dp^.  The  time-domain  Internal  magnetic  field 
resulting  from  an  incident  EMP  signal  has  been  calculated  for  some 
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representative  cases.  Its  time  dependence  can  be  understood  on  the  basis 
of  the  "low-pass"  character  of  both  the  composite  shell  material  and  the 
cylindrical  geometry  Itself. 


The  principal  difficulty  In  treating  screened  boron-epoxy  composite 
laminates  lies  In  finding  an  appropriate  sheet-lmmlttance  description  for 
the  screen.  We  have  developed  such  a description  appropriate  for  the  case 
where  the  mesh  dimensions  are  electrically  small,  and  which  takes  the 
presence  of  the  composite  panel  Into  account.  Our  result  constitutes  an 
extension  of  the  classic  work  of  Kontorovlch.  Using  this  method  of 
characterizing  the  mesh  we  have  addressed  the  problems  of  plane-wave 
transmission  through  a planar  screened  boron-epoxy  composite  layer  and 
penetration  Into  the  interior  of  a cylindrical  shell  of  screened  composite 
material.  We  have  found  that  the  composite  layer  itself,  except  for  its 
effect  on  the  sheet  impedance  of  the  screen.  Is  relatively  unimportant 
insofar  as  the  shielding  behavior  of  the  screened  laminate  is  concerned. 
Thus  one  needs  to  consider  only  the  sheet  impedance,  properly  modified  to 
take  into  account  the  presence  of  the  composite  layer,  in  the  analysis  of 
shielding  problems. 

The  principal  difference  between  graphite  and  screened  boron-epoxy 
composites  from  a signal-transmission  standpoint  is  that  the  graphite 
composite  tends  to  act  as  a low-pass  or  integrating  filter,  while  the 
screened  boron-epoxy  composite  behaves  as  a high-pass  or  differentiating 
filter.  Thus  the  time-domain  field  transmitted  through  a planar  layer  of 
screened  composite  contains  a contribution  proportional  to  the  derivative 


of  the  input  waveform,  as  a result  of  the  inductive  component  of  the  screen 
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response  dominated  by  this  derivative  term,  while  the  later-time  response 
is  influenced  most  strongly  by  the  wire  conductivity  in  the  mesh. 

The  cylindrical  screened-composite  shell  is  of  particular  Interest 
from  an  analytical  standpoint  because  the  Integrating  effect  of  the 
cylindrical  geometry  tends  to  cancel  the  differentiating  effect  of  the 
wire-mesh  screen.  In  fact,  for  an  ideal  screen  with  perfectly  conducting 
wires  the  Internal  magnetic  field  waveform  is  identical,  except  for  a 
reduction  in  amplitude,  to  the  incident  waveshape.  For  the  screen  para- 
meters which  have  been  considered  in  this  report,  it  has  been  shown  that 
the  energy  shielding  ratio  principal  break  frequency  tends  to  be  lower 
than  that  of  a comparable  graphite  composite  shell,  especially  when  the 
screen  wires  are  highly  conductive.  Thus  the  screened  boron-epoxy  com- 
posite shell  can  be  a much  better  shield  than  a graphite  composite  shell  as 
a consequence  of  its  larger  equivalent  conductivity-thickness  product. 
However,  we  have  also  noted  that  the  equivalent  conductivity-thickness 
product  alone  is  not  a completely  valid  descriptor  of  screened  composite 
panels;  it  is  only  useful  at  low  frequencies,  or,  equivalently,  at  late 
times.  At  early  times  the  mesh  size  is  the  critical  parameter,  since  it 
determines  the  Inductance  of  the  sheet.  In  fact,  we  have  presented 
numerical  data  which  indicate  that  overall,  a screen  with  a smaller 
equivalent  conductivity-thickness  product  can  be  a better  shield  than 
another  with  a larger  product,  at  least  with  respect  to  the  incident-field 
waveforms  considered  here.  This  result  suggests  the  possibility  that  one 
might  even  design  a screen,  whose  wires  have  a given  conductivity  and 
which  possesses  a given  amount  of  metal  per  unit  area,  for  optimum  shielding 
performance  against  a given  incident  waveshape. 
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2.  CONCLUSION 

It  is  clear  that  while  advanced  composite  laminates  are  not  yet  as 
effective  as,  say,  aluminum  panels  in  shielding  against  unwanted  penetra- 
tion of  electromagnetic  fields,  they  are  nonetheless  good  enough  to  be  used 
as  aircraft  or  missile  skins  without  placing  an  unreasonable  burden  on 
the  internal  shielding  (conduit,  caHe  braids,  etc.)*  This  fact,  coupled 
with  the  strength  and  weight  advantages  of  advanced  composites,  guarantees 
the  increased  future  use  of  these  materials,  and  points  up  the  importance 
of  continued  study  both  of  their  intrinsic  electrical  properties  and  of 
the  fundamental  electromagnetic  boundary-value  problems  which  involve  them. 

Some  suggested  problems  for  future  study  include  the  following: 

(a)  Integral-equation  formulation  of  scattering  problems  involving  ^ 
hollow  composite  shells,  using  the  boundary  connection  operator 

(b)  Singularity  Expansion  Method  (SEM)  solutions  of  transient 
scattering/penetration  problems  involving  follow  composite 
shells  of  separable  or  non-separable  shape 

(c)  Study  of  effects  of  Imperfect  panel-to-panel  joints  and  joints 
between  metallic  and  composite  panels 

(d)  Field  penetration  through  apertures  in  composite  panels 

(e)  Continued  study  of  bonded  and  unbonded  wire  meshes,  including 
the  effects  of  linear  and  nonlinear  junction  impedance  and 
dielectric  environments 

(f)  Canonical  problems  Involving  the  propagation  of  surface  waves 
on  wire  meshes,  including  interactions  with  nearby  conductors; 
composite  waveguides. 

These  are  only  a few  of  the  many  problems  which  come  to  mind.  Advanced 
composites  provide  a solution  to  many  problems  because  of  their  great 
strength  and  light  weight;  but  to  the  electromagnetic  theorist  and  the 
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engineer  concerned  with  electromagnetic  shielding,  they  are  a source  of 
many  future  challenges.  It  Is  sincerely  hoped  that  this  report  will 
contribute  to  an  Improved  understanding  of  the  electromagnetic  shielding 
behavior  cl  these  materials  and  excite  the  Interest  of  other  investigators 
who  will  continue  to  study  them. 
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